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Q-l Abstract. Pangenesis is the mechanism for jointly producing the visible and dark mat- 

ter asymmetries via Affleck-Dine dynamics in a baryon-symmetric universe. The baryon- 
symmetric feature means that the dark asymmetry cancels the visible baryon asymmetry 
and thus enforces a tight relationship between the visible and dark matter number densities. 
The purpose of this paper is to analyse the general dynamics of this scenario in more detail 
and to construct specific models. After reviewing the simple symmetry structure that under- 
pins all baryon-symmetric models, we turn to a detailed analysis of the required Affleck-Dine 
dynamics. Both gravity-mediated and gauge-mediated supersymmetry breaking are consid- 
ered, with the messenger scale left arbitrary in the latter, and the viable regions of parameter 
space are determined. In the gauge-mediated case where gravitinos are light and stable, the 
regime where they constitute a small fraction of the dark matter density is identified. We 
discuss the formation of Q-balls, and delineate various regimes in the parameter space of the 
Affleck-Dine potential with respect to their stability or lifetime and their decay modes. We 
outline the regions in which Q-ball formation and decay is consistent with successful pangen- 
esis. Examples of viable dark sectors are presented, and constraints are derived from big bang 
nucleosynthesis, large scale structure formation and the Bullet cluster. Collider signatures 
and implications for direct dark matter detection experiments are briefly discussed. The 
following would constitute evidence for pangenesis: supersymmetry, GeV-scale dark matter 
mass(es) and a Z 1 boson with a significant invisible width into the dark sector. 
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1 Introduction 

The present-epoch matter content of our universe points to the existence of new particles 
and interactions beyond those described in the Standard Model (SM). Precision cosmological 
measurements have established that 4.6% of our universe consists of visible matter (VM), 
and 23% consists of dark matter (DM). 

The known properties of the visible matter assert that its relic abundance today can be 
understood only if we hypothesise yet-unknown processes that created an asymmetry between 
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the visible baryons and antibaryons in the early universe. The dynamical generation of such 
an asymmetry requires processes which violate baryon number, C and CP conservation, and 
occur out of equilibrium at earlier times. The survival of the baryonic asymmetry today 
implies that any baryon-violating processes must have been ineffective thereafter. 

On the other hand, the properties of dark matter are poorly understood, and there are 
many possible mechanisms for the production of the DM relic density. Well-studied DM 
scenarios rely on extensions of the SM which are motivated by other unsolved problems in 
particle physics. Examples of DM candidates include: weakly interacting massive particles 
(WIMPs), expected to arise in TeV extensions of the SM (in particular the LSP in super- 
symmetric models); the axion, introduced to solve the strong CP problem; sterile neutrinos, 
hypothesised to explain the neutrino masses; mirror DM in a parity-symmetric universe. 

It is intriguing, however, that despite the possibly disparate origin of the visible and 
dark matter relic densities, and thus the different parameters determining their values, these 
densities differ only by a factor of a few. The apparent coincidence between the VM and DM 
is reinforced by the (as yet unverified) observations of the DM direct-detection experiments 
DAMA, CoGeNT and CRESST, which favour DM masses of a few GeV. This mass scale, 
which is very similar to that of the nucleons, may be considered together with the similar 
visible and dark total mass abundances to suggest that the number densities of the visible 
and dark relic species are related. It is then important to consider how such a relation can 
arise dynamically. 

A tight connection between the VM and the DM densities is established dynamically 
if the visible and the dark sectors are both charged under a common particle number which 
remains always conserved, and under which the two sectors develop compensating asym- 
metries. The common and always conserved symmetry has to be a generalisation of the 
ordinary baryon number, and we may thus refer to this paradigm as a "baryon-symmetric 
universe" . The separation of the baryonic and antibaryonic charge into the visible and the 
dark sectors, respectively, amounts to an asymmetry generated in an appropriately-defined 
particle number orthogonal to the generalised baryon number [1] . We describe the symmetry 
structure and symmetry breaking pattern that gives rise to a baryon-symmetric universe in 
Sec. 2. It should be already clear though, that the separation of baryons and antibaryons 
into two sectors requires non-equilibrium dynamics. 

While establishing such a connection between the VM and DM densities is motivated 
phenomenologically, it is also desirable that this occurs within extensions of the SM which 
solve other (fundamental) particle-physics problems. Reference [1] proposed that a baryon- 
symmetric universe can arise in super symmetric models (SUSY) from Affleck-Dine (AD) 
dynamics [2, 3]. This mechanism was termed "pangenesis", signifying the simultaneous 
generation of all matter. Examples of pangenesis were outlined in [1] and subsequently in [4], 
although explicit dark sectors were not constructed. 

In this paper, we discuss the mechanism of pangenesis in more depth, and construct 
specific models. We consider both gravity- or Planck-mediated SUSY breaking (PMSB) and 
gauge- mediated SUSY breaking (GMSB). In the latter case, we analyse the Affleck-Dine dy- 
namics for arbitrary messenger mass, and present semi-analytical results which are of general 
relevance for AD baryogenesis in SUSY models, not just for those embodying pangenesis. We 
estimate the asymmetry generated, which in pangenesis quantifies the separation of baryonic 
and antibaryonic charge, and discuss the formation and decay of Q-balls, the non-topological 
solitons which can arise via the AD mechanism. We produce plots which exhibit the feasibil- 
ity of successful AD pangenesis (or baryogenesis) . We also construct a class of simple dark 
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sectors and describe the cascade of the dark baryonic charge down to the lightest dark sector 
particles. We discuss the nature of DM (which here can be 'atomic'), and various pertinent 
constraints. 

The possibility of a baryon-symmetric universe was first discussed by Dodelson and 
Widrow [5-7], and subsequently developed in various models which employed different dy- 
namics for the separation of the baryonic from the antibaryonic charge. Besides the AD 
mechanism, these are: out-of-equilibrium decays of heavy particles [8-16], with such sce- 
narios termed "hylogenesis" [14]; the QCD phase transition [17]; asymmetric freeze-out [18]; 
asymmetric freeze-in [19]; a time-dependent background of a light scalar field [20], and bubble 
nucleation due to a first-order phase transition [21]. 

Related to the baryon-symmetric proposal, but distinctly different both in symmetry 
structure and observational signatures, is the scenario in which a baryonic asymmetry is 
originally generated in either the visible or the dark sector, and shared at a later time between 
the two sectors via chemical equilibrium. This possibility has been explored in Refs. [22-75]. 
Moreover, VM and DM may share the baryonic asymmetry in comparable amounts, without 
the need of a separate dark sector, if the baryonic asymmetry arises via the AD mechanism 
and is partially stored in the form of stable Q-balls comprising the DM of the universe [76]. 
In this class of models, there is no symmetry, common to both the visible and the dark sector, 
which remains conserved throughout the cosmological history of the universe. However, both 
in this and in the baryon-symmetric scenario, the relic abundance of DM is due to a particle 
asymmetry, and both scenarios are thus often referred to as 'asymmetric DM (ADM)' models. 

This most generic property of baryon-symmetric models - that there is an always con- 
served particle-number symmetry - implies the possible existence of an associated gauge 
boson. The latter can be discovered in colliders and provides a channel for DM direct detec- 
tion. This probe is absent in models with no unbroken baryonic symmetry. 

Finally, we note that it is possible to relate the visible and the dark matter relic abun- 
dances without relying on a dark asymmetry. The possibility that the AD mechanism is 
responsible for both the visible baryonic asymmetry and a symmetric DM abundance has 
been discussed in Refs. [77-91]. These models consider non-thermal DM production from 
late Q-ball decays. Possible connections between the WIMP miracle and baryogenesis have 
been explored in Refs. [92-94]. 

In the next section we describe the symmetry structure and symmetry breaking pattern 
that give rise to a baryon-symmetric universe. We discuss their realisation in supersymmetric 
models via the Affleck-Dine mechanism - the process of pangenesis - in Sec. 3. In Sec. 4, 
we expound on the technical aspects of the AD dynamics. We construct a specific dark- 
sector model in Sec. 5, and discuss various cosmological constraints in Sec. 6. We discuss 
observational signatures of pangenesis in Sec. 7. 

2 How to build a baryon-symmetric universe 
2.1 Symmetry structure 

The relic particle content of our universe is determined largely by the symmetries of the 
low-energy effective interactions, as described by the SM in the visible sector. 

Most of the present VM mass density is carried by protons and heavier nuclei, whose 
stability is understood in terms of a global U(l) symmetry of the SM, generated by the baryon 
number operator of the visible sector, which we shall call B v . The charge B v is additively 
conserved in perturbation theory, though violated by non-perturbative effects associated with 
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anomalies. These non-perturbative effects, which are very weak for temperatures below that 
of the electroweak phase transition, violate the global By 1 , but conserve a combination of the 
baryon and lepton symmetries of the SM, (B — L) v . The conservation of the global {B — L) v , 
and of B v at low temperatures, ensures that an existing excess of baryons over antibaryons in 
the visible sector can be sustained, as indeed observed in today's low-energy universe. Most 
ways of understanding the origin of the asymmetry between visible baryons and antibaryons 
require that the (B — L) v symmetry was violated at some stage in the early universe, by 
as-yet unknown processes which subsequently became completely ineffective. 2 

If the DM relic abundance is also attributed to a particle-antiparticle asymmetry, then 
at least some of the dark-sector species have to be charged under a particle number which 
is additively conserved at low energies. We shall refer to this symmetry as the 'dark baryon 
number' B d . In analogy to the visible sector, a dark baryonic asymmetry must have originated 
from processes which once violated B d under the appropriate conditions for creation of a net 
B d charge. Subsequently, these processes must have become unimportant, restoring B d as a 
good symmetry of the low-energy interactions, and maintaining the existing dark baryonic 
asymmetry. 

The visible and the dark baryonic asymmetries will be related if the processes which gen- 
erated them violated the {B — L) v and B d symmetries while respecting a linear combination 
of the two. Quite generally, we may consider the two linear combinations 3 



If the 'generalised baryon number' B — L remains conserved throughout the cosmological 
evolution of the universe, while X, C and CP violating interactions occur out of equilibrium, 
a net X charge will be created, amounting to related (B — L) v and B d charge asymmetries, 



where r][Q] = Yli1i[ n (fi) ~ n (fi)]/ s i s the charge-to-entropy ratio of charge Q. The visible 
and the dark baryonic asymmetries compensate each other under the generalised baryon 
number B — L, and thus constitute a 'baryon-symmetric universe'. The X asymmetry quan- 
tifies the separation of baryonic and antibaryonic charge into the visible and dark sectors, 
respectively. The restoration of the X symmetry in the low-energy late universe, together 
with the conservation of the B — L, ensure that the visible and the dark baryonic asymmetries 
are preserved separately. 

The strict conservation of B — L is guaranteed if it is a gauge symmetry. A gauged B — L 
has multiple functionality for pangenesis. It precludes the possibility of the AD mechanism 
being operative along flat directions of the scalar potential with non-vanishing B — L, thus 
ensuring the validity of Eq. (2.2). It facilitates the decay of the heavier of the visible and 
dark sector LSPs, thus reducing a heavy relic that is unwanted in our scenario. Moreover, 

x Note that they respect a discrete subgroup of it, and thus do not endanger proton stability. 

2 For completeness we note that the visible-sector baryon asymmetry can also be explained through elec- 
troweak baryogenesis, which uses the non-perturbative (B — L) v -conserving but £? v -violating 'sphaleron' pro- 
cesses mentioned earlier. 

3 Under the definitions of Eq. (2.1), visible-sector fields are those with B — L — X, and dark-sector fields 
are those with B — L — —X. There may also be connector-type fields with linearly independent B — L, X. 



B - L = (B - L) v - B d 
X = (B-L) v + B d . 




r ] [(B-L) v ]=rj[B d ] 



V [X] 



(2.2) 
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the associated gauge boson offers an important probe into the dark sector, in both collider 
and DM direct-detection experiments. We discuss these points in subsequent sections. 

However, a gauged B — L symmetry has to be broken to be consistent with the non- 
observation of a long-range force. We shall require that the breaking of the local B — L 
leaves a global symmetry unbroken, under which all of the degrees of freedom appearing in 
the low-energy effective interactions transform according to their gauged B — L charges. We 
discuss the spontaneous breaking of B — L in the next subsection. 

2.2 Gauged and spontaneously broken B — L 

We wish to spontaneously break the gauged XJ(1)b-l symmetry without spoiling the symme- 
try structure discussed in the preceding subsection. In particular, we must avoid introducing 
deleterious side-effects such as too-rapid proton and/or DM decay, and asymmetry erasure. 
The first step to ensuring that these goals are met is to construct the gauged U(1)b-l as 
a diagonal subgroup of U(l)( B „ i ) v xU(l)B d xU(l)s, where this last symmetry acts only to 
transform the Higgs field(s) responsible for the breaking of the gauged U(l)#-L- In other 
words, we seek to identify the gauged B — L generator as 

(B - L) gaugcd = (B - L) v - B d + S = (B - L) global + S, (2.3) 

with S, (B — L) v and B d [hence also (B — L) v — B^ = (B — L) g i bai] being purely global 
symmetry generators when acting through any linear combination orthogonal to the one 
above. 

We now explain this in more detail. Suppose, for definiteness, that the visible sector 
is simply the minimal supersymmetric standard model (MSSM). Then, as per the previous 
discussion, (B — L) v is just the usual generator for global B — L transformations of the 
quark and lepton superfields; in particular, it transforms neither the dark sector fields nor 
the Higgs sector responsible for breaking the gauged U(l) b-l- Similarly, B<± generates an 
analogous global phase symmetry in the dark sector, and does not act in either of the other 
two sectors. The third generator S transforms only the Higgs supermultiplets responsible for 
the breaking of the local U(1)b-l; it does not transform any visible or dark sector fields. A 
simple example of such a sector is exhibited through the superpotential terms 

WD PT(aa - v 2 B _ L ), (2.4) 

where a and a constitute a vectorlike pair of chiral superfields, and T is trivial under all 
symmetries. The U(l)g global symmetry transforms a and a with equal and opposite phases, 
with 5-charges q a and —q a , respectively. These superpotential terms contribute to the scalar 
potential the term 

V D /3 2 \aa - v 2 B _ L \ 2 (2.5) 

which induces nonzero VEVs, (era) = v B _ L , for the scalar components of the superfields 
a, a. We denote the scalar component of a superfield by the same symbol as the superfield, 
allowing context to resolve the ambiguity. 

This VEV spontaneously breaks the global symmetry U(l),g and the gauge symmetry 
U(l)(£-L) gaugcd without breaking the global symmetries generated by (B — L) v and B&. The 
potential Nambu-Goldstone boson from spontaneous U(l)s- breaking is eaten by the B — L 
gauge field, so there is no massless scalar boson in the physical spectrum. It is worth noting 
that requiring the MSSM sector to respect (B — L) ga ugcd; and hence automatically (B — L) v , 
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ensures that standard global baryon number and, separately, global lepton number emerge 
as accidental symmetries at the renormalisable level. 

To ensure that U(l)g exists as an independent symmetry, we must choose the charge 
q a judiciously. For example, we should not choose q a = 1, for that would allow the superpo- 
tential terms LH u a and n c a without violating overall B — L conservation. 4 (The superfield 
a would then essentially contain a left-handed antineutrino.) While we could impose inde- 
pendent U(l)s invariance even in this case, it is much more desirable to have this symmetry 
arise accidentally. At the renormalisable level, this is ensured in the visible sector provided 
q a 7^ ±1/2, ±1, ±2 so that neither a nor a can couple to the monomials LH U , (n c ) 2 and n c . 

We wish to extend this desirable feature to non-renormalisable effective operators as 
well. The discussion below will focus on operators connecting the visible sector to the o- 
sector, but an extension to dark sector fields will be required once that sector is specified. 

Returning to the MSSM-gauge-invariant monomials LH U and n c , we now need to con- 
sider higher order terms given generically by LH u a k and LH u a k and similarly for n c . All 
such terms will be automatically absent if q a ^ ±l/k for any positive integer k. At cubic 
order, the MSSM monomials u c d c d c , LLe c and Qd c L all have (B — L) v = —1, so the same 
requirement forbids their couplings to positive-integer powers of a and a also. At higher 
order, in addition to other MSSM monomials carrying {B — L) v = ±1, there are monomials 
carrying (B — L) v = ±2 and —3 [95]. Some of these monomials are listed in Table 2. By mul- 
tiplication, one can combine these 'basis monomials' into a countable infinity of terms which 
must have (B — L) v charges of the form ±n =t 2m — 3p, where n, m and p are non-negative 
integers. The (B — L) v charge of the general term O n)rri)V a k is ±n ± 2m — 3p + q a k, where 
m ^p is any visible-sector monomial. For any rational charge q a , there exist monomials that 
can couple to a positive integer power of a or a in a [B — L)-invariant way. Unless we allow 
the unorthodox possibility of making q a irrational, it is thus not possible for U(l)s to be an 
accidental symmetry to all orders. 5 Nevertheless, it is obvious that one can choose a rational 
q a to forestall the explicit breaking of U(1)(£-l) v xU(1)s to U(l)(£-L) v +s to arbitrarily high 
order in terms of mass dimension. Thus, we conclude that the spontaneous breaking of the 
gauged U(l)#-L poses no in-principle problem for proton decay and baryon asymmetry era- 
sure. By a straightforward extension, we will also be able to ensure sufficient DM stability 
and preservation of the dark-sector asymmetry. 

Neutrino masses do not alter the above picture if neutrinos are Dirac particles. This 
includes the possibility of Dirac seesaw, which could explain the smallness of neutrino masses. 
If neutrinos are Majorana, then limits on the Majorana masses -and thus the B — L breaking 
scale- apply, in order to avoid wash-out of the visible and dark baryonic asymmetries. These 
limits depend on the seesaw scheme. 

3 Pangenesis via Affleck-Dine 

In the Affleck-Dine mechanism [2, 3], the coherent oscillations of a scalar field generate a 
charge under a U(l) symmetry, when this symmetry is conserved by the low-energy cou- 
plings but explicitly broken by high-energy interactions. The formation of a condensate - a 

4 We denote the chiral superfields of the MSSM by the standard notation Q, L, u c , d c , e c and n c for 
the quark doublet, lepton doublet, up antiquark, down antiquark, charged antilepton and neutral antilepton, 
respectively. 

5 This conclusion has been reached assuming the visible sector is described by the MSSM. If the visible- 
sector gauge symmetry is extended, then the form of the gauge-invariant effective operators will be more 
constrained. 
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homogeneous field configuration that possesses a time-dependant phase corresponding to a 
net U(l) charge - relies on the scalar field acquiring a large vacuum expectation value (VEV) 
in the early universe. A large VEV amplifies small U(l)- and CP-violating terms, making 
them important at early times. If such terms are present, and if the scalar field possesses 
a large VEV but is not confined to a static configuration, i.e. it is free to spontaneously 
relax towards a lower energy state, the conditions for creation of an asymmetry [96] are sat- 
isfied, and a U(l) charge is generated. Subsequently, the expansion of the universe causes 
the amplitude of the oscillating field to redshift, thus effectively restoring the U(l) symmetry 
and maintaining the generated charge. This charge may be transferred to lighter degrees of 
freedom transforming under the same symmetry if the condensate decays or thermalises via 
U(l)-conserving couplings. 

The above dynamics can be naturally realised in supersymmetric theories with low- 
energy global U(l) symmetries. Because of the high degree of symmetry, the scalar potential 
of SUSY models generically possesses many flat directions (FDs) at the renormalizable limit. 
The FDs are typically lifted by non-renormalizable operators, which may violate the U(l) 
symmetry that characterises the FD at low field values. The effective potential along a 
FD contains also SUSY-breaking contributions induced by the hidden sector, the thermal 
bath and the vacuum energy of the universe. The latter contribution dominates during 
inflation, and in the epoch of inflaton oscillations before reheating is completed, during 
which H ^> m s , where H is the Hubble parameter and m s stands for the soft mass scale 
induced by the SUSY-breaking hidden sector. The large vacuum energy of the universe and 
the absence of (unsuppressed) quartic terms along a FD can naturally drive the FD field to a 
large VEV at early times. However, as H decreases, a vanishing VEV becomes energetically 
more favourable. When the FD field begins to relax towards a zero VEV, the U(l) and CP 
violating terms induce a time dependent phase for the field, which now thus carries a non-zero 
U(l) charge. If the U(l) symmetry of the FD is a symmetry of the full low-energy effective 
theory, the condensate decay and thermalisation processes will respect it, and transfer this 
charge to the lightest particles carrying the same conserved number. 

This is the case in the MSSM, where most of the FDs of the scalar potential carry baryon 
or lepton number. The AD baryo/lepto-genesis in the MSSM generates a (B — L) v charge, 
which is carried initially by a squark or slepton condensate, and is subsequently transferred 
to the fermionic degrees of freedom by the {B — L) v -conserving interactions of the MSSM. 
The FDs of the MSSM have been listed in Ref. [95]. 

Pangenesis can be realised in extensions of the MSSM along FDs with vanishing gener- 
alised B — L charge and non-zero X charge: 

D B _ L = 0t Tfl _ L0 = O ( 3. 1} 
D x = $T X <\> + , (3.2) 

where T B _ L and T x are the B — L and X generators, respectively. The field (ft parametrises 
the FD, and is in general a linear combination of visible, dark and connector sector fields. A 
sum over these components of <j) is implied in Eqs. (3.1) and (3.2). The X charge of the FD 
field is 

q XA = D x /\ ( f>\ 2 . (3.3) 

Equation (3.1) is warranted along FDs if the generalised B—L is gauged. In this case, the 
corresponding L>-terms in the scalar potential lift the directions which carry non-vanishing 
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B — L, thus asserting that the AD mechanism cannot generate a net B — L charge. 6 While 
a gauged B — L is most attractive for ensuring that it is always conserved, it is possible to 
justify a zero B — L otherwise. In particular, a minimal Kahler potential for field directions 
with non-zero B — L induces large positive mass-squared contribution along these directions 
during the vacuum-energy dominated era. This ensures that such field directions do not 
acquire large VEVs, and thus renders the AD mechanism inoperative along them. 

Equation (3.1) can be satisfied along field directions with either T B L cp = 0, or T B _ L cp 7^ 
and D B L = 0. In the former case, the FD consists of connector-sector field(s) which carry 
both visible and dark baryonic charge, with (B — L) v = B^ ^ 0. An example of this was 
presented in [1]. The latter case corresponds to field directions along which B — L is spon- 
taneously broken while its expectation value remains zero. Such directions are necessarily 
made up by multiple fields, which may be purely visible and purely dark sector fields (car- 
rying only visible or dark baryonic charge), or may also involve connector-sector fields with 
various charge assignments. Examples of this case were presented in [4]. 7 

In the next section, we present explicitly the FD potential which gives rise to the AD 
dynamics, in the PMSB and the GMSB scenarios. We present extensions of the MSSM in 
which these dynamics operate under the conditions of Eqs. (3.1) and (3.2) in Sec. 5. 

4 The Affleck-Dine Dynamics 

The directions of the scalar potential which are flat in the renormalisable and supersymmetric 
regime are typically lifted by non-renormalisable and SUSY-breaking interactions. 
We consider non-renormalizable terms in the superpotential of the type 

which lift a FD parameterised by the scalar degree of freedom <I>. It will often be convenient 
to reparametrise the scalar component, denoted also by <I>, in terms of amplitude and phase, 
as per 

$ ee ^=e ie . (4.2) 

The dimension d of the operator in Eq. (4.1) is determined by the particle content and the 
charge assignments of the model under consideration. This term violates explicitly the U(l) 



6 We note that, for field_yEVs below the (B — L)-breaking scale M B _ L , these D-terms are suppressed by 
a factor ^1 + M 2 _ L /m 2 ^ , where m s is the soft scale (see e.g. [97]). To ensure that the AD field does not 
deviate from a path of vanishing B — L during the important stages of its dynamics, we shall assume that 
M B L is smaller than the VEV <j> oa c of the AD field at the onset of oscillations (respectively the VEV at 
Q-ball formation if applicable). Since the (B — L)-gauge boson can be as light as 500 GeV for a B _ L ~ 
(see Sec. 7.1), this can always be satisfied. 

7 We note though, that the two examples of Ref. [4] are likely not viable: They employ dimension-3 and 
dimension-4 monomials to lift the FDs and generate the asymmetry. It is known, however, that dimension-4 
monomials along a single flat direction do not produce enough asymmetry due to thermal corrections [98-100] . 
For a dimension-3 monomial, this problem can only become more severe. In addition, a dimension-3 monomial 
means that the direction considered is not .F-flat already at the renormalisable level, and an unnaturally small 
coupling is required in order to generate sufficient asymmetry. This is an unnecessary complication arising 
from the specific choice of gauge-invariant monomial. Indeed, the scalar potential of supersymmetric theories 
typically possesses a large number of directions which are both D-flat and F-flat, at least at the renormalisable 
regime. This is in fact one of the basic reasons for which the Affleck-Dine mechanism is implemented in 
supersymmetric models only: that the quartic couplings in the scalar potential vanish identically along flat 
directions, at the supersymmetric regime, due to the symmetries of the theory. 
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symmetry characterising the FD at low field values. For successful pangenesis, this symmetry 
has to be U(l)x- The scale M* at which U(l)x is violated may be as high as the Planck 
scale, and as low as the requirement of non-erasure of the baryonic asymmetries allows. 8 We 
examine this issue in Sec. 5. The superpotential coupling of Eq. (4.1) contributes to the 
potential along the FD the term 



VsuSY 



Qd-t 

— . (4.3) 



Despite arising from an A-violating coupling in the superpotential, this term respects the 
U(l) symmetry of the FD. It dominates and stabilises the potential at large field values. 

SUSY breaking also contributes to the potential along FDs. There are three sources 
which are relevant for the AD dynamics: the vacuum energy of the universe, the thermal 
bath and the hidden sector which dominates SUSY breaking in the late universe. 

During inflation and in the epoch of inflaton oscillations, the universe is dominated by 
the energy stored in the inflaton field. The coupling of the inflaton to the AD field can induce 
a large effective mass for the latter. The mass term depends on the Kahler potential: while 
a minimal coupling gives rise to a positive mass-squared term, non-renormalisable couplings 
in the Kahler potential can generate a negative mass-squared contribution for the AD field, 

V H = -cH 2 \<S>\ 2 , (4.4) 

where c ~ 1. A detailed study of the Hubble-induced terms for various types of the Kahler 
potential can be found in [103]. A Hubble-induced negative mass-squared contribution en- 
sures that the FD field acquires a large VEV at early times, which is essential for the efficacy 
of the AD mechanism. The VEV of the field is determined by the balance between the terms 
of Eqs. (4.3) and (4.4): 

rh • / // \ • 

*M TT ) ■ (4-5) 



y/2 \M*, 

The SUSY breaking in the hidden sector also induces a soft mass term along the FD, 
whose precise form depends on the way SUSY breaking is communicated. The mass-squared 
contribution from the hidden sector has to be positive at low energies (small VEVs) to ensure 
the restoration of the symmetry 

V s = m s 2 |$| 2 + ... , (4.6) 

where the dots stand for corrections which depend on the SUSY-breaking mediation mech- 
anism, and which become important with increasing field VEV. In the following sections, 
we consider both the PMSB and the GMSB scenarios, and specify the form of the term of 
Eq. (4.6). 

In addition, the dilute thermal bath of the universe during the epoch of inflaton oscilla- 
tions generates a thermal potential for the AD field [98, 99]. The details of thermalisation of 
the inflaton decay products depend largely on the model of inflation, and can be affected by 
large VEVs developed along flat directions of the scalar potential [104, 105]. Here we assume 
that the inflaton couples generically to both the visible and the dark sector, and model the 
temperature of the inflaton decay products before the completion of reheating by [106] 

T « (HTlM P ) 1/4 , (4.7) 



Note that this scale can be effectively larger than the Planck scale when the operator arises within a 
Froggatt-Nielsen framework (see e.g. [101, 102]). 
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where Tr is the reheating temperature. This gives rise to two contributions: (i) The rela- 
tivistic degrees of freedom which couple directly to the AD field induce a thermal mass term 
for the latter: 

V T Dc k y 2 k T 2 \<S>\ 2 , for y k <P<T. (4.8) 

Here c k are constants of order 1, and y k stand for the Yukawa or gauge coupling constants 
of the AD field to these particles, defined such that y k 4> is the effective mass they pick up 
when the AD field has a non-zero VEV (p. A summation over all species for which y k <p < T is 
implied in Eq. (4.8). The degrees of freedom for which this condition does not hold, obtain a 
large mass and become non-relativistic. The suppression of their contribution to the thermal 
mass of the AD field at y k (f>/T > 1 is exponential, (ii) If some of the heavy degrees of freedom 
for which yicfi > T, also couple to other relativistic particles, they can be integrated out, and 
generate effective couplings of the AD field to the thermal bath. The large (^-dependant mass 
of these mediators affects the running of their couplings to the light particles, and results in 
a potential which is logarithmic in the amplitude of the AD field [99] 

VxD Si if^ r4ln (^-) ' for ^> T > ( 4 - 9 ) 

where again q = 0(1) and summation over all degrees of freedom for which y\§ > T is 
implied. The l/167r 2 -factor reflects the loop suppression of this contribution. 9 
The effective potential felt by the AD field at low energies is 

V eS = V s + V T . (4.10) 

The interplay between the hidden-sector /thermal and the Hubble-induced mass terms de- 
termines the spontaneous relaxation of the AD field towards its late-time ground state, 
as governed by its equation of motion $ + 3H& = —dV/d&*. When H becomes suffi- 
ciently small, the hidden-sector /thermal mass term drives the evolution of the field because 
V^ s ((j)) > V' H (4>), and the field starts oscillating around the origin. This occurs at 

Hox ^ s (Ykm v \ (4 .n) 



The amplitude of the field at the onset of oscillations, ijl osc , is given by Eq. (4.5) with 
H = H osc . It is during this non- equilibrium process of relaxation towards the lowest energy 
state that the asymmetry is generated. 

For this to occur, it is of course necessary that U(l)-breaking terms operate along 
the FD. Such terms are generated as a result of SUSY breaking 10 and the U(l)-violating 
superpotential interaction of Eq. (4.1): 

Va = A — — j— o + h.c. (4.12) 



9 We do not use the interpolation between Eqs. (4.8) and (4.9) which was introduced in Ref. [107], 
Vt ~ T 4 ln(l + y 2 |$| 2 /T 2 ), since it overestimates the contribution of Eq. (4.9) to the effective mass-squared 
of the condensate by 16n /yf. This factor is significant even for very large couplings, yi ~ 1. This fact is 
important in determining the onset of oscillations of the AD field and thus the generated asymmetry (see 
Eqs. (4.11) and (4.19)). 

10 In AD scenarios in which multiple flat directions are excited, U(l) violation arises also in the terms 
derived directly from the superpotential. This enhances the asymmetry produced [108-110]. This feature was 
employed in the model of pangenesis constructed in [1]. 
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In the above, the parameter A is again determined by the mediation mechanism of SUSY 
breaking to the FD fields. There may also be a Hubble-induced contribution, 5A ~ H. 
However, such a contribution is absent in .F-term inflation models, in which the interaction 
leading to it is proportional to the inflaton field and thus averages to zero in the epoch of 
inflaton oscillations [111]. We shall omit this contribution in our analysis. 

We ultimately want to estimate the A-charge density generated in the condensate via 
the AD mechanism, 

n x = iq x - , (4.13) 

where q x is the X charge of the FD field. The equations of motion in an FRW background 
give 

n x+ 3Hn x = igx ^*-^ . (4.14) 

Only the A-violating terms of the FD potential contribute to the right-hand side of Eq. (4.14). 
In the setup we consider here, X violation arises only from the A-terms of Eq. (4.12), thus 
Eq. (4.14) becomes 

\A\ U/V2) d , , s 

n x + 3Hn x = 2q x ' ' v ^_ 3 ; sin (9 A + d ■ 9) , (4.15) 

where we parametrise A = \A\ exp(i0A), and <j) and 9 are defined in Eq. (4.2). The above 
equation reflects two essential features of the AD mechanism: (i) the large VEV of the AD 
field amplifies the effect of the otherwise suppressed U(l)-violating interactions; (ii) the CP 
violation necessary for the generation of an asymmetry arises from a misalignment between 
the phase of an initial field configuration and the phase of the U(l)-violating terms. 

The CP violation is the final element we need to identify in this set-up. At early times, 
the field is dynamically confined at a large VEV by the Hubble-induced negative mass term, 
which generates a curvature of order H along the radial direction, around the minimum. The 
A-term generates d minima and maxima along the angular direction in the internal phase 
space of the field, with radial curvature of order A. If A <C H, the angular curvature of 
the potential is much smaller that the radial curvature, and the field picks a random initial 
phase. 11 When the radial oscillations begin, the mismatch between this initial phase and the 
phase of the A-term kicks the field in the angular direction, thus generating a time-dependant 
phase and a net charge for the condensate. 

The charge is thus generated at the onset of oscillations, within a Hubble time, while 
CP violation is effective. It can be estimated from Eq. (4.15), by ignoring the instantaneous 
expansion of the universe: 5n x ~ h x 5t with 5t ~ H~ , which gives 



-(t osc ) ps 2q x sin<5 ^ J^° S 3 C ^^ w q x sin S \A\^ ac , (4.16) 



where in the second equality we used Eq. (4.5). The phase S represents the effective CP 
violation, and we shall typically assume sin 5 ~ 1. 



n If A does not contain a Hubble-induced contribution, this is typically the situation. However, if the 
hidden-sector potential is very flat, oscillations may begin very late, and A ~ ff sc can be realised, even in 
the absence of Hubble-induced ^-terms (see Sec. 4.2). 
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The asymmetry generation occurs in the epoch of inflaton oscillations, during which 
there is rapid entropy production. The X-charge-to-entropy ratio should then be calculated 
at reheating, after which the entropy per comoving volume remains constant. The A-charge 
density decreases as n x oc a -3 . During inflaton oscillations, the Hubble rate scales like 
H oc a~ 3 / 2 . This remains valid even after the inflaton has decayed if most of the energy is 
stored in non-relativistic particles. This can happen if large VEVs of scalar fields after infla- 
tion induce large masses for the inflaton decay products [104, 105]. As these VEVs redshift, 
various degrees of freedom become light, and contribute to the relativistic energy density 
of the universe and to the various thermalisation processes. Assuming that the transition 
from a matter-dominated-like universe to the radiation-dominated era signals the complete 
thermalisation of all degrees of freedom (reheating), the A-charge density at reheating is 
estimated to be 12 

H 2 

n x {t R ) =n x (t osc )-f . (4.17) 

-"osc 

The entropy density at reheating is s(£r) = AH R Mp /Tr, thus the X-charge-to-entropy ratio 
is 

n x (t R ) _ n x (t osc )T R 
Vx ~ s(t R ) ~ 4iJ 2 sc M| ' (4 ' i8j 

Using the expression on the right-hand side of Eq. (4.16), this gives 

_ q x sm5 T R \A\(t> 2 QSC ^ g x sm5 \A\T R ( M* \ ^ , . 

~ 4 Mp H 2 SC ~ 2 M 2 \H os J ' l4 - i9j 

where in the second equality we made use of Eq. (4.5). We will use both forms of Eq. (4.19) 
in our analysis below. 

An important quantity for the late-time evolution of the condensate is the ellipticity of 
the field orbit in its internal phase space: 

9 n v /q v „ \A\ 

e = ~ = - JLl Tr =sin<5^--L . 4.20 

LO LO(j) z H osc 

The quantity uj, defined in Eq. (4.11), is the maximum rotational velocity that the AD 
field can develop at the onset of oscillations, 9 ^ uj. Maximum 9 occurs if the ^4-terms 
which provide the U(l) violation are dominant, A = H osc , and the CP-violation is maximal, 
sin 5 ~ 1. 

4.1 Gravity-mediated SUSY breaking: asymmetry generation 

We shall now specify the form of Eq. (4.6) in the PMSB scenario, combine all the contributions 
to the FD potential discussed above, estimate the asymmetry generated, and identify the 
parameter space for successful pangenesis. 



12 There may be a period of radiation domination between this transition (at Hubble rate _ff m _ r ) and 
complete thermalisation (at Ha). For H osc > ff m _ r > Hr, Eq. (4.17) should in this case be replaced 
by n x (t R )/n x (t osc ) = [H m - r /H osc ] 2 [Hr/ 'H m - r ] 3/2 = [H R /H osc ] 2 [H m - r /H R } 1/2 . This means that the X- 
charge-to-entropy ratio becomes enhanced by a factor of [H m - T / Hr] 1 ^ 2 compared to Eq. (4.19). Considering 
the various upper bounds on the reheating temperature (determined by Hr) from cosmological considerations, 
this opens additional parameter space for successful asymmetry generation via the AD mechanism. In the 
rest of the analysis, we shall only focus on the -ff m _ r = Hr case, as described in the text. 
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The potential along a FD in the PMSB scenario can be described by [112] 



V AD =m 2 s \<f>\ 2 
where 



1 + iHn 



|$|2 

m 2 



V T -cH 2 \$\ 2 + [A T -*+h.c. + 



d-l 



M, 



d-3 



(4.21) 



A « m s , (4.22) 



and the soft scale is of order the gravitino mass in PMSB, m s ~ ^3/2 • The first term in 
Eq. (4.21) is the hidden-sector SUSY breaking mass term. The correction proportional to K 
arises radiatively at 1-loop, and depends on the FD [112, 113]. Typically \K\ ~ 0.01 — 1. The 
sign of K can be important for the late-time evolution of the AD condensate, as we discuss 
in Sec. 4.3. The thermal potential Vt is given by Eqs. (4.8) and (4.9). 

The AD field starts oscillating around the origin of the potential when, according to 
Eq. (4.11), the Hubble expansion rate becomes 

~ 4 

tt2 ^ 2 1 2 rr2 1 c lVl 1 osc / A 00 \ 

H osc ~ m s + c k y k T osc + 16?r2 ^ 2 , (4.23) 

where summation over all degrees of freedom for which yk<p < T and yi<f> > T is implied. 
Equations (4.5), (4.7) and (4.23) have to be solved together to obtain the time at which 
oscillations start. Typically, the thermal contributions to the potential force oscillations 
to start earlier, and thus suppress the asymmetry generation, as per Eq. (4.19). This, in 
turn, means that a larger reheating temperature and/or a larger M* are required in order 
to obtain a sufficiently large X-charge-to-entropy ratio. The precise effect of the thermal 
corrections on the onset of oscillations and the asymmetry generation depends non-trivially 
on M*, Tr and the various coupling constants y^, and more detailed analyses can be found in 
Refs. [101, 114]. Here we intend to only roughly sketch the parameter space which produces 
the correct asymmetry. We thus pick a sample value for a coupling, y ~ 10~ 2 , calculate the 
onset of oscillations numerically, and present the combination of M* , Tr values which yield 
■q x ~ 10 -9 in Fig. I. 13 

The reheating temperature is constrained from the production and decay of gravitinos 
in the early universe. Gravitinos are produced from scatterings of thermal particles in the 
primordial plasma. Even if they are unstable, which is typically the case in PMSB scenarios, 
their lifetime is extremely long: a gravitino lighter than 10 TeV decays at time scales larger 
than 1 sec, after BBN has started. The successful predictions of BBN can then be retained 
only if the abundance of gravitinos is sufficiently small. This yields an upper limit on the 
reheating temperature. A recent detailed analysis was presented in Ref. [117] which finds 
that for an unstable gravitino of about 500 GeV mass, the reheating temperature has to 
satisfy T R < 10 6 GeV. 

Given this constraint, Eqs. (4.19) and (4.23) imply that for d = 4 the asymmetry 
generated is not sufficient, even for M* as high as Mp (this is true even in the limit y — > 0). 



13 For the purpose of illustration, in all of the graphs we use m s = 500 GeV. The equations of course 
retain their explicit dependence on the soft scale. While LHC has already placed stringent bounds on the 
masses of the superpartners, it is still possible for some of the sparticles to be at or below 500 GeV [115, 116]. 
In addition, pangenesis occurs along flat directions which involve both MSSM and dark-sector scalars. The 
masses of the latter are not constrained, and may be lower than in the visible sector. A lower soft scale in 
the dark sector may in fact be desirable: it can be generated dynamically, and may be responsible for the 
GeV-scale DM mass (see discussion in Sec. 5.4.2). The soft scale relevant for mixed visible-dark flat directions 
is an average of the visible and the dark soft scales, and can be thus comfortably below 500 GeV. 



- 13 - 



10 



,19 




M, > M P 



T R (GeV) 



Figure 1. The Tr — M* values which yield the correct asymmetry in PMSB. The different curves 
correspond to various values of the dimension d of the monomial which lifts the flat direction. From 
top to bottom: d = 5 (blue), d = 6 (red), d — 7 (yellow), d = 8 (green). The dashing of the 
lines denotes the excluded region from the requirement that Q-balls decay before the LSP freezes 
out, T Q > 10 GeV, which sets T R , > 6 • 10 3 GeV (see Eq. (4.60)). This bound does not apply if 
Q-balls do not form, as is the case for flat directions with K > 0, and even for K < it may 
be relaxed (see text for discussion). The vertical dotted line at Tr = 10 6 GeV denotes the upper 
limit on the reheating temperature in PMSB scenarios, from considering the production and decay of 
gravitinos [117]. The horizontal dotted line is drawn at the reduced Planck mass. For this plot, wc 
used r/ x = 10- 9 , q x = 1, sin S = 1, y= 10~ 2 , m s = 500 GeV and K = -0.01. 

However, for d ^ 5, it is possible to get rj x ~ 10 -9 for M* < Mp while respecting the upper 
limit on the reheating temperature, as shown in Fig. 1. 

4.2 Gauge-mediated SUSY breaking: asymmetry generation 
4.2.1 Effective potential and onset of oscillations 

We now specify the hidden-sector SUSY-breaking mass term of Eq. (4.6), in the GMSB 
scenario, present the complete potential, and follow the introductory analysis in order to 
estimate the asymmetry generation and investigate the parameter space. 

The potential along a flat direction due to GMSB from the hidden sector has been 
calculated in the limit of small and large field VEVs in Ref. [118]. Here, we adopt the 
interpolation introduced in Ref. [88]. Including the terms already discussed, we consider the 
following potential 



where m s is the soft mass scale and M m is the messenger scale. 

14 In Eq. (4.24), we have omitted the contribution from PMSB which becomes important only at large field 
amplitudes, when ml /2 (f> 2 > mlM^. Using Eq. (4.32), we find that this occurs at <j> > KT 3 Af P ~ 10 15 GeV, 
independently of both the gravitino mass and the soft scale. However, the analysis of Sec. 4.2.3 shows that 
the oscillations of the AD field in the GMSB scenario always start at lower field VEVs. It is thus justified to 
omit the PMSB contribution in the subsequent analysis. 




(4.24) 
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The first term in Eq. (4.24) resembles a quadratic potential, m^ 2 /2, in the <p <S M m 
limit, and it flattens out for <j) > M m . This is because in this regime, the gauge fields acquire 
large masses ~ g<j) > M m and the transmission of SUSY breaking is suppressed. As in the 
PMSB scenario, the correction proportional to K arises from 1-loop corrections, and depends 
on the FD. 

The j4-terms in GMSB are generated at 2-loop and are thus suppressed with respect to 
the soft scale. Following [88], we model the parameter A in Eq. (4.12) by 



a m s 

A = m 3/2 + — TT/2 ■ 425 

' W 1 



Here the first term is the contribution from PMSB, and the second term includes the loop- 
suppression factor, a Q ~ 0.01, and takes into account that GMSB is suppressed at large field 
VEVs 4> > M m . 

For the potential of Eq. (4.24), the oscillations around the origin begin according to 
Eq. (4.11) at 

/ H \ 2 In ( 1 + -4r?h) 

(4 26) 



)1K 



1 + 



where we omitted corrections which are proportional to K and logarithmic in cf) osc /^/2M m , 
and also the thermal corrections. The latter can be unimportant in the regimes allowed by 
other constraints for reasonable values of the FD couplings to other fields. We discuss this 
further in Sec. 4.2.3. At the onset of oscillations the field VEV is determined by the minimum 
of the potential prior to oscillations, Eq. (4.5), which for later convenience we rewrite as 

^ *b , (4.27) 



\/2M m V m . 

where we defined 

'mMi 



M d-2 



l 

d-2 



(4.28) 



Equations (4.26) and (4.27) yield the amplitude of the field, 4> osc , and the Hubble 
parameter, H osc , when oscillations begin. The dimensionless parameter b defines different 
parametric regimes, which we list in Sec. 4.2.3. In each case, we estimate the generated 
asymmetry from Eq. (4.19). The asymmetry depends on the reheating temperature Tr, 
which in the GMSB scenario is constrained from the requirement that gravitinos do not 
overdose the universe. We now turn to this issue. 

4.2.2 Constraints from gravitinos 

The gravitino couples to the other species with strength inversely proportional to its mass. 
In GMSB, gravitinos are expected to be light and stable, thus scatterings and decays of 
thermalised supersymmetric particles in the primordial plasma can produce a significant 
gravitino relic abundance [118, 119]. 

For m 3 /2 ^ 100 keV, scattering processes produce gravitinos more efficiently than de- 
cays, and result in a gravitino relic density 

^-<^)(i^) 2 (i^)' 
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where m g \ is the gluino mass. We shall require that gravitinos are a subdominant contribution 
to the dark matter of the universe, i.e. f2 3 / 2 < 0.1 f^M- This yields an upper limit on Tr, 
which becomes less severe if the gravitino is heavy. 

For smaller masses, 100 keV > ^3/2 ^ 1 keV, decays of supersymmetric particles pro- 
duce gravitinos very efficiently. It is then necessary to ensure that squarks and sleptons are 
never very abundant, which implies that the reheating temperature has to be lower than 
their mass. To be conservative, we will require Tr < 100 GeV in this regime. 

For even lower masses, 1 keV > 7713/2 ^ 100 eV, if the reheating temperature is too high, 
the gravitinos thermalise, decouple while relativistic, and overdose the universe. We therefore 
require again Tr < 100 GeV to ensure that the relevant scattering processes are never in 
(relativistic) equilibrium, and gravitinos do not thermalise despite their large coupling to the 
MSSM particles. If m 3 / 2 < 100 eV, the gravitino relic abundance after their freeze-out is 
subdominant, and there is no limit on the reheating temperature. 

While a larger gravitino mass reduces the gravitino relic abundance (which is unwanted 
for our scenario), thus allowing a higher reheating temperature without overproducing gravi- 
tinos, it also decreases the decay rate of the NLSP, since the gravitino couples more weakly. If 
the NLSP was thermalized in the early universe and the gravitino is too heavy, the late decay 
of the NLSP freeze-out abundance can spoil BBN. We therefore again require Tr < 100 GeV 
in this regime. The value rh of the gravitino mass, for which this bound becomes relevant, 
depends on the nature and the mass of the NLSP. It is more stringent for bino than for stau, 
and quite relaxed for sneutrino NLSP [117]: 

1 GeV, for NLSP = b 
fh={ 10 GeV, for NLSP = f (4.30) 
100 GeV, for NLSP = v . 

In summary, we consider the following constraints on the reheating temperature: 

' 100 GeV, for m 3 / 2 > fh 



Tr<{ 



2 



m 3/2 „ ... (4 31) 



100 GeV, for 100 eV < m 3/2 < 100 keV 

no limit, for ?n 3 / 2 % 100 eV , 



where the constraint of the first line arises from BBN considerations, while the rest ensure 
underabundant gravitino DM. 

Finally, we note that the LSP, instead of being the gravitino, could also reside in the 
dark sector. We outline a model where this is possible in Sec. 5.4.2. The gravitino then 
decays to this lighter particle. This would ease the bound on the reheating temperature 
from the requirement that the LSP is a subdominant contribution to the dark matter. In 
addition, the NLSPs in the visible sector decay to the dark sector via the (B — L)-gaugino, 
resulting typically in a shorter NLSP lifetime than in the case of gravitino LSP. This would 
ease constraints from NLSP decays during BBN. A more detailed study of the constraints 
on the reheating temperature in this case is, however, beyond the scope of this work. 
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Figure 2. The field amplitude over the messenger scale, <f> osc / \/2M m , and the Hubble parameter over 
the soft mass, H osc /m s , at the onset of oscillations, plotted as a function of the dimensionless variable 

b = (m s M?- 3 /M?- 2 ) 1/(d ~ 2) for the GMSB potential of Eq. (4.24). The various curves correspond to 
different values of the dimension d of the monomial which lifts the flat direction: d = 5 (solid, blue), 
d = 6 (dashed, purple), d — 7 (dot-dashed, yellow), d = 8 (dotted, green). The behaviour of the 
curves is described by Eqs. (4.33) and (4.39). 



4.2.3 Parameter space analysis 

In what follows we estimate the charge-to-entropy ratio generated via the AD mechanism in 
the GMSB scenario, according to what was described in Sec. 4.2.1. We apply the constraints 
discussed in Sec. 4.2.2, and identify the available parameter space. 

The charge-to-entropy ratio, tj x , depends on four parameters: the reheating temperature 
Tr, the messenger mass M m , the scale M* at which the U(l)x violation occurs, and the 
dimension d of the polynomial which lifts the FD. We take m s ~ 500 GeV for definiteness. 
The requirement that r\ x ~ 10 -9 reduces the free parameters to three. This means that 
for a fixed value of d, the required reheating temperature is a function of M* and M m , i.e. 
Tr = J^(M*,M m ). The reheating temperature is independently constrained, as summarised 
in Eq. (4.31), and the constraints depend on the gravitino mass. In GMSB, the gravitino 
mass is related to the messenger scale by 

4vr V rn s ) \ 

It is then convenient, by means of Eq. (4.32), to express the reheating temperature that is 
required for the generation of sufficient asymmetry, in terms of M* and the gravitino mass, 
Tr = ^(M*, m 3 / 2 ). This makes it possible to map the reheating-temperature constraints 
on the 777.3/2 — M* or the M m — M* plane, for various values of d. We do so in the graphs 
of Fig. 3, where we also sketch various regimes associated with the production and decay of 
Q-balls, which we discuss in Sec. 4.3. 

This requires to first determine the onset of oscillations of the AD field. Equations (4.26) 
and (4.27) can be solved for the dimensionless variables (j) osc /^/2M m and H osc /m s , in terms 
of the dimensionless parameter b, defined in Eq. (4.28). We present the numerical solution 
in Fig. 2, and we shall now describe the semi-analytical behaviour. For completeness, we 
also discuss why the thermal corrections of Eqs. (4.8) and (4.9), which were neglected in 
Eq. (4.26), are not relevant in the parameter region allowed by gravitino constraints. 

We discern the following cases: 



(4.32) 
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For b < 1, the solution to Eqs. (4.26) and (4.27) is well approximated by 



V2M n 



1 , 



(4.33) 



and A ~ a m s + ni 3 / 2 , where a Q ~ 0.01 is the loop-suppression factor of the A-term in 
GMSB. Since ^> osc < V2M m , the oscillations begin when the field is in the quadratic part 
of the potential. The A-term is subdominant, and the ellipticity of the condensate is 



~ sin — 



sin 5 



a Q + 



< 1 . 



As per Eq. (4.19), the charge-to-entropy ratio is 



q x sin 5 ( A \ 2 / M£ Tr 



m s M 2 



(4.34) 



(4.35) 



osc with Tosc;. Let us first show that </> sc 



Thermal corrections: We need to compare 
thus the thermal corrections are of the type of Eq. (4.9) - rather than Eq. (4.8) - for rea- 
sonable values of the couplings y. In what follows, we will adopt the value m s « 500 GeV, 
to make the estimates more definite. 

At H = H osc ~ m s , the temperature is T osc ~ (rngT^Mp) 1 / 4 , where Tr is determined from 
Eq. (4.35) with the requirement that rj x = 10~ 9 . Setting q x = sin 5 = 1, this gives 

T osc w 10 21 GeV 2 /0o 



(4.36) 



The significance of the thermal corrections is maximal for minimal (j) osc . We will now deter- 
mine the minimum value of <^ OS c allowed by the constraints on the reheating temperature. 
Combining Eqs. (4.31) and (4.35), we find 



V2 bM m > I 



10 12 GeV 



1 GeV 



3/2 



1/2 



4 • 10 14 GeV, 



if < fh 

if m 3 / 2 > rh 



(4.37) 



where rh is defined in Eq. (4.30). Note that because we are in the b < 1 regime, Eq. (4.37) 
also implies that M m > 10 12 GeV or m 3/2 > 0.5 GeV. 

Comparing Eqs. (4.36) and (4.37), we see that OSC 3> T osc . We shall thus consider the 
thermal potential of Eq. (4.9). The thermal correction to H^ sc arising from Eq. (4.9) is of 
the order [(y 2 /4ir) T 2 sc /<^ osc ] 2 . Let us now compare this to the value of H osc ~ m s given in 
Eq. (4.33). Using the above, we get 



y 2 T 2 C /0 OSC 



4tt 



< 



3/2 



10- 6 y 2 , 



if m 
if m 



3/2 
3/2 



< m 



> m . 



(4.38) 



Thus we find that, if m 3 /2 < fh, the thermal contribution to H osc is subdominant at least 
f° r V ^ 0.1(1 GeV/m 3 / 2 ) 3 / 4 - If ^3/2 > on t ne other hand, the thermal correction to 
-f/osc is obviously irrelevant. 
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For 1 < b < [rn s /m^i2) d ~ 2 , the solution to Eqs. (4.26) and (4.27) is well approximated by 



V2M n 



b d -i . 



(4.39) 



Since </> osc > \/2M m , the oscillations begin when the field is in the flat part of the potential. 
Using the above solution, the parameter A of Eq. (4.25) becomes 



A ~ to 3 /2 + a m s b d - 1 
The ellipticity of the condensate then is 



e b>1 =sm6— 



A 



sin 5 



(4.40) 



(4.41) 



and remains less than unity since b^ 1 < m s /m 3 / 2 in this regime. Using Eq. (4.19), the 
X-charge-to-entropy ratio is 



q x sin (5 / A \ <d-2) / M^T R 



11 K 



— } b d-l 



m s Mp 



(4.42) 



Thermal corrections: We shall again compare 4> osc with T osc . At H = H osc ~ m s (M m /(j) osc ), 
the temperature T osc « (.HoscT^Mp) 1 / 4 is 

T osc « 10 21 GeV 2 (M 3 /^L) 1/4 ~ 4 • 10 30 GeV 2 (m| /2 /^ sc ) ^ , (4-43) 

where we solved Eq. (4.42) for Tr, requiring as before that t/ x = 10~ 9 . Using the con- 
straints of Eq. (4.31) on Tr, we find that 4> osc must satisfy 



10 12 GeV, 



V2b^ M m > 



10 14 GeV 



m 3 / 2 \V3 



1 GeV 



if m 3 /2 < ^ 
, if m 3 / 2 > m . 



(4.44) 



Equations (4.43) and (4.44) can be combined to show that 4>osc/T osc 3> 1 for the range of 
gravitino masses of interest. Thus, we shall again consider the thermal correction to H osc 
arising from Eq. (4.9), and compare it to the value of H osc given in Eq. (4.39). Using the 
above, we get 



y 2 ^osc/ 'Pose ^ ^2 



Air K 



(JSC 



6 ■ 10 12 GeV 



\ 7/2 / ^3/2 y/ 2 

J U GeV/ 



(4.45) 



Applying Eq. (4.44), we find 



y 2 T 2 C /^ < 



47T H osc ~ 



60092 (rHk) 1 ' 2 ' if 

1 GeV^ 2/3 



m 3 / 2 < m 



6 • 10~ 5 y 2 



m 3/2 ) 



(4.46) 



, if m 3 / 2 > m. 



If 777,3/2 < rh, the thermal contribution is subdominant at least for y < 0.01, and may only 
enhance H osc by a factor of a few for 0.1 < y < 1. For m 3 / 2 > ?7i, on the other hand, the 
thermal correction to H osc is completely negligible for the allowed values of Tr. 
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d-1 



• For b > (m s /m3/ 2 ) d ~ 2 3> 1, the oscillations begin in a region of the potential where 
the contribution from PMSB is relevant: The corresponding term in the potential is 
Vpmsb ~ m 3/2l^ ) l 2 ■^ c l - Comparing with the first term in Eq. (4.24), we see 

that PMSB effects are important at the onset of oscillations if m^^&osc ~ m s^m- From 
Eq. (4.39) (which is valid up to the overlapping region between the two regimes), we find 

d—i 2) 

that this happens for b > (m s /m 3 / 2 ) d ~ 2 or, equivalently, (m 3 y 2 M t 3 ) > 10 _3 Mp. 

Requiring that M* < Mp and m 3 / 2 < m s, this regime is only possible for d > 8. We will 
restrict ourselves to d < 8 when presenting the allowed regions of parameter space and 
have checked that this regime is excluded by constraints from Q-ball decay for d = 8 (see 
Fig. 3). We shall therefore not consider this case further. 

Using the above analysis of the dynamics which generates the asymmetry, we numeri- 
cally translate the bounds of Eq. (4.31) into exclusion regions on the m 3 / 2 — M* plane, for 
d = 5 — 8 and present the results in Fig. 3. For each parameter pair in the allowed regions, the 
X-charge-to-entropy ratio 7] x ~ 10~ 9 can be obtained for a value of the reheating tempera- 
ture which satisfies the bounds in Eq. (4.31). Additional constraints for successful pangenesis 
arise from considering the late-time evolution of the AD condensate, which we now discuss. 

4.3 Formation and decay of Q-balls 

The Affleck-Dine condensate is unstable with respect to spatial perturbations [76, 78, 107, 
112, 113, 120-124]. If the potential of the AD field grows more slowly than 2 , the AD 
condensate can fragment into bound states of the scalar fields participating in the flat di- 
rection, known as Q-balls [125]. Q-balls are non-topological solitonic configurations of fields 
which transform under a global U(l) symmetry. They carry a non-zero charge under this 
symmetry, and they are stable against decay into quanta of the same field. If they form, they 
may affect the cascade of the charge generated in the condensate down to the lightest visible 
and dark-sector baryons. 

Q-balls can decay into gauge-invariant combinations of particles which carry the same 
charge. They may decay into scalar degrees of freedom that are linearly independent from 
the Q-ball field, and/or fermions. A given Q-ball decay channel is kinematically allowed only 
if the mass-to-charge ratio of the Q-ball is larger than that of the decay products. The decay 
into fermions proceeds only from the surface of the Q-ball, while the decay into scalars can 
occur from the interior. The latter is thus typically enhanced with respect to the former. 
We parametrise this enhancement of the decay rate into scalars by a factor /q which can 
be as large as Jq < 10 3 , depending on the flat direction [78]. For decays to fermions, we 
correspondingly set /q = 1. The decay rate of Q-balls has an upper limit of [126] 

where loq is the rotational speed of the Q-ball field configuration, <3?q = [<pQ(r) / y/2\ exp(iwQt), 
in the internal phase space, Rq is the radius and Q is the charge of the Q-ball. Note that 
the charge is measured in units of 4> quanta and is related to the X charge by 

Qx = q x Q (4.48) 

where q x is the X charge assignment of the FD, defined in Eq. (3.3). The upper limit on the 
decay rate is reached if y<f)n /ojq 3> 1 , where y is the coupling of the Q-ball field to the fields 
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Figure 3. Various regimes on the m 3 / 2 — M* or M TO — M* plane for asymmetry generation via the 
Affleck-Dine mechanism in GMSB models, for arbitrary messenger scale. Here, m 3 / 2 is the gravitino 
mass, M m is the messenger scale, and M* is the scale of the operator which lifts the flat direction. 
The scenario of pangenesis is implemented as described in this paper, within the non-shaded areas 
for different values of the reheating temperature. 

The meaning of the shaded areas is as follows (cf. Eqs. (4.3f) and (4.85) - (4.88)): 

(i) In the blue area, gravitino DM is dominant or overabundant. 

(ii) In the orange areas, the NLSPs decay late and may spoil BBN. The upper bound on m 3 / 2 depends 
on the nature of the NLSP; it is at I GeV for bino, 10 GeV for stau and 100 GeV for sneutrino NLSP. 
Above the orange line, the NLSPs do not thermalise while relativistic, their abundance is always 
suppressed, and their decay does not affect BBN. 

(iii) In the red area, the Q-balls formed after the fragmentation of the AD condensate are stable. 

(iv) In the yellow area, Q-balls decay at low temperatures Tq < 10 MeV, during or after BBN. 

(v) In the green area, the gravitinos produced in the decay chain Q ball — > NLSP — > gravitino may 
become overabundant (we note that this is a stringent limit, and could be significantly relaxed). 

(vi) In the gray area, the scale of the operator that lifts the FD and produces the asymmetry exceeds 
the reduced Planck mass, M* > Mp. 

To guide the eye, we have drawn the lines 6=1 and b*^ = 70, where b = (m s M*~ 3 /M^ -2 ) d 2 \ 
For b < 1 (dashed line on the right), the oscillations of the AD field start in the quadratic part 

d-2 

of the potential, and the Q-balls formed are gravity-type. For b d - x > 70 (dashed line on the left), 
Q-ball decay does not produce NLSPs (except at the late stages of evaporation). The four graphs 
correspond to different dimensionality of the operator which lifts the flat direction and is responsible 
for the U(l)x violation: d = 5 upper left, d = 6 upper right, d — 7 bottom left, d = 8 bottom right. 
For these graphs, we have used m B = 500 GeV, t] x = 10~ 9 , q x = 1, sin<5 = 1 and /q = 1. 
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it decays into. The decay rate is typically saturated for Q-balls formed via fragmentation 
of an AD condensate and which couple to visible and dark sector fields (see e.g. [90]). The 
quantities Rq, ojq and cf)Q depend on the parameters of the low-energy AD potential and 
the charge Q. However, the charge Q itself is determined from the fragmentation of the AD 
condensate, and thus depends on the dynamics of the AD mechanism at higher energies. 

The fragmentation of an AD condensate and the formation of Q-balls has been studied 
semi-analytically and numerically. In what follows, we use these studies to delineate various 
regimes in the parameter space of the AD potential, with respect to the formation, stability 
or lifetime, and decay modes of the Q-balls. We depict these regimes in Fig. 3. 

For the purposes of the pangenesis scenario described in this paper, we shall focus on 
the regions where if Q-balls form, they are unstable and decay sufficiently early, releasing all 
of their charge to the thermal plasma without spoiling any of the predictions of the standard 
cosmology. To be conservative, we shall require that the Q-balls formed after Affleck-Dine 
pangenesis decay before BBN, the latest at temperature 

T Q > 10 MeV . (4.49) 

We outline this region in the graphs of Fig. 3, using the properties of Q-balls described in the 
following subsections. Additional stronger requirements may arise, depending on the type of 
Q-ball, which we shall now discuss. 

4.3.1 Q-balls in gravity- mediated SUSY breaking 

After the onset of oscillations, the dominant part of the FD potential for the case of PMSB 
in Eq. (4.21) is 

1 + K In 



$1 



with m s ~ "1,3/2 ~ 500 GeV. The formation of Q-balls depends critically on the sign of the 
parameter K: If K > 0, the potential grows faster than quadratic, and no Q-balls are formed. 
If K < 0, the potential grows more slowly than quadratic, and Q-ball solutions exist. 

The correction proportional to K arises at 1-loop, and typically \K\ ~ 0.01 — 1. Gaugino 
loops contribute 5K < 0, while Yukawa couplings give 5K > [113]. In the MSSM, K > 
arises typically for purely leptonic FDs (since they do not couple to the gluino), and for 
baryonic FDs with large third-generation contribution (as long as the gluino is lighter than 
the stops). For baryonic FDs of first and second generation fields, K < 0. 

If K < 0, most of the charge density carried initially by the condensate gets trapped in 
Q-balls. The average charge of these Q-balls can be estimated from numerical simulations 
as [124] 

Q « 2-10- 2 f^V « 4-10- 2 (MlX^ , (4.50) 

where we used Eq. (4.5) and the fact that H osc ~ m s in PMSB, and assumed ellipticity e ~ 1. 
For the potential in PMSB, the mass, radius and rotational speed of a Q-ball with charge Q 
are [78, 112] 

M Q = m s Q , (4.51) 

= W^ s ' (452) 

ujq = m s . (4.53) 
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Using the above and Eq. (4.47), the temperature at the time of Q-ball decay is 



For the charge of Eq. (4.50), we then find 
' f Q m s M P ^ 



5 x /4 ff Q m s M P 
\K\Q 



1/2 



(4.54) 



Tr 



\K\ 



m s 



d-3 
d-2 



Jq_ q x sin 6 T R \ 



1/2 



2\K\ 



M P J 



(4.55) 



where in the second equality we used Eq. (4.19) to express M* in terms of the reheating 
temperature Tr (and the X-charge-to-entropy ratio r] x ). 

The mass-to-charge ratio of Q-balls in PMSB, Mq/Q ~ m s , allows them to decay into 
fermions and generically also into scalars (unless the FD field consists only of the LSP). 
If Q-balls decay after the LSP freezes out, they will introduce a non-thermal LSP density, 
which has to be subdominant for successful pangenesis. If this density is sufficient to bring 
the LSPs in equilibrium, annihilations will reduce it to 



n 



Q 

LSP 



(Tq) 



H 

(av) 



(4.56) 



T=Tc 



The abundance determined by Eq. (4.56) is the upper limit on the LSP relic density from Q- 
ball decays. We emphasise that the actual abundance may be smaller. It is useful to compare 
^lspCTq) with the LSP relic density from thermal freeze-out (which has to be a subdominant 
component of DM in models of pangenesis, as we discuss in Sec. 5). At temperature T = Tq, 
the latter is 



n 



Lp(Tq) ^ n[ SP (T f ) ^ c ^ 



T 3 
1 f 



T 3 



H 

(av) 



(4.57) 



T=T f 



where Tf is the freeze-out temperature of the LSP. Combining Eqs. (4.56) and (4.57), we get 



n 



Q 

LSP 



71 



LSP 



(Tq) 



n 
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LSP 



71 



/ 

LSP 



(Tq) 



Tj_ (av) f 

t q {™)q 



Tf Tf 



(4.58) 



for s- or p-wave annihilation of the LSPs. If f^ L SP ^ ^dm, it is possible to accommodate 
Q-ball decay after the LSP freeze-out while maintaining subdominant LSP abundance. For 
simplicity and in order to avoid any extra assumptions, however, we shall require Tq > Tf. 
Typically, Tf « ?tt,lsp/20, and we thus impose 



T Q >T f ~ 10 GeV 



(4.59) 



Given Eq. (4.55), the constraint of Eq. (4.59) poses an upper limit on M* or, equivalently, 
a lower limit on Tr: 



1 



T R > 10 4 GeV — 



\K\ 



fQ V0.01 



500 GeVV 



77l s 



T 



10 GeV 



(4.60) 



We denote this limit in Fig. 1 (for the conservative choice of /q ~ 1). We note, however, 
that this is a worst-case-scenario limit if Q-balls indeed form (K < 0), which does not apply 
if Q-balls do not form (K > 0). 
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4.3.2 Q-balls in gauge-mediated SUSY breaking 

After the onset of oscillations, the dominant part of the FD potential for the case of GMSB 
in Eq. (4.24) is 



Fad = m 2 M 2 In 2 1 + 



l + K\n 



Following the analysis of Sec. 4.2.3, we discern two cases: b < 1 and 1 < b < (m s /m 3 / 2 ) 
where b = (m s Mf" 3 /M^ 2 ) 1,{d ~ 2) , as defined in Eq. (4.28). 



d-l 
d-2 



If b < 1, the oscillations of the AD field start in the quadratic part of the potential, at 
</> osc /\/2 ~ bM m , and the situation resembles the PMSB case. Q-balls form only if K < 0. If 
they form, most of the charge density of the condensate is carried by Q-balls of charge [124] 

' rn 3/2 n2 ^500 GeV\ 
1 GeV/ 



Q 



b<l 



10" 



10 17 b 2 



(4.61) 



where the prefactor is suppressed with respect to Eq. (4.50) due to the lower ellipticity of the 
condensate, e < 1, in the GMSB case. The other properties of the Q-balls are as described 
in the previous subsection, Eqs. (4.51) to (4.54). Substituting Eq. (4.61) into (4.54), we find 
the temperature at the time of Q-ball decay 



Q,b<l 



250 GeV 



1 (I GeV\ 



m 



3/2 



1 



0.01 

wr 



500 GeV 



(4.62) 



Because of their large mass-to-charge ratio, Mq/Q ~ m s , the decay of these Q-balls 
produces NLSPs, which subsequently decay into gravitinos and imparity even fermions. We 
must ensure that the gravitinos remain a subdominant component of DM. A calculation 
analogous to that of Eqs. (4.56) to (4.58) gives the maximum gravitino relic density as 
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3/2 
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3/2 Tf {av)f ^ m 3/2 
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Tf T 2 
To^T 2 



(4.63) 



where ^nlsp ^ s * ne w °uld-be NLSP thermal-relic density if NLSPs did not decay, and the 
two cases in brackets are for s- and p-wave annihilation, respectively. Typically Tf ~ m s /20. 
Assuming f^LSP ~ ^dm and requiring figL < ^dm, we get a rough constraint on Tq. For 
the case of s-wave annihilation, this gives 



Tq > m 3/2 /20 . 

Given Eq. (4.62), this constraint can be written in terms of b and m.3/2 as 



m 3/2 \ 
1 GeV/ 



,3 A ^°- 01 
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\K\ 



500 GeV 



(4.64) 



(4.65) 



Since & < 1, this constitutes a rather mild upper limit on 7713/2, even for /q = 1. 



d ~ 2 , the oscillations of the AD field start in the nearly- 

' d-2 J 

~ })d-i > 1. Numerical simulations show that most of 



For the case 1 < b < (m s /m 3 /2} 
flat part of the potential, (p osc /^/2M n 
the charge density is trapped inside Q-balls of average charge [107] 

f 6 2 ^ 2 



(4.66) 
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where (3 depends on the ellipticity of the condensate: (3 ~ 6 • 10 4 for e ~ 1 and /3 ~ 6 ■ 10 
for e < 0.1. Thus 
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(4.67) 



1 GeW \ m s J 

Note that the two formulas for the Q-ball charges in the GMSB case, Eqs. (4.61) and (4.67), 
agree reasonably well in the overlapping region of validity b ~ 1, given that they are based on 
two different numerical simulations (Refs. [124] and [107], respectively). The mass, radius, 
rotational speed of the field and field amplitude at the center of Q-balls which form in the 
nearly-flat part of the potential are [76, 120, 121, 127, 128] 
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Q 1 / 4 
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4> Q = {m.M^ 2 Q X ' A . 
The mass-to-charge ratio of this type of Q-ball thus decreases with charge: 
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(4.72) 



500 GeV 



For sufficiently large Q, the mass-to-charge ratio is smaller than the lightest gauge-invariant 
combination of quanta carrying unit charge. In this case, Q-balls are energetically disallowed 
to decay, and in many scenarios found in the literature they constitute (part of) the DM 
of the universe (see e.g. [76, 122, 129-133]). Q-ball dark matter carrying (visible) baryonic 
charge is viable as long as the Q-balls were formed along FDs lifted by visible-baryon-number 
violating operators [128]. 

For the pangenesis scenario described here (see Sec. 5), we require that the Q-balls 
formed can decay, so that all of the asymmetry generated in the AD condensate is released 
into the plasma, and cascades down to the lightest visible and dark sector particles (we will 
explore the possibility of antibaryonic Q-ball dark matter in the framework of pangenesis 
elsewhere). This sets an upper bound on the charge of the Q-balls formed. A more stringent 
bound arises from requiring Q-balls to decay sufficiently early. In this regime, the temperature 
at the time of Q-ball decay is 
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(4.73) 



The minimum temperature at which the Q-balls decay without changing the predictions of 
standard cosmology depends on their decay products. These are determined by the mass-to- 
charge ratio Mq/Q. There are the following possibilities: 



- 25 - 



Large Q-balls with mass-to-charge ratio 

— < ^ f i ~3GeV (4.74 
Qx Q 22fQxj 

are stable. Here ^/ m / ^ s the sum of masses of the lightest gauge-singlet combination of 
fields which carry a net X charge ^2fQ Xf - In pangenesis, this combination involves both 
visible and dark sector fields with equal X charges, q Xvis = q x dark , such that a (B — L)- 
singlet combination is formed. The lightest dark-sector combination of fields carrying 
1^0 charge, and zero net charge under any dark-sector forces (except B — L), is the DM 
state itself. As we shall see in Sec. 5, in order to get the right relic DM mass density, 
the mass of the DM state (which may be composite) has to be m DM ~ 5 GeV x q x DM (cf. 
Eq. (5.35), for asymmetry release in the thermal bath below the electroweak phase tran- 
sition). The visible-sector mass contribution to the right-hand side of Eq. (4.74) depends 
on the FD along which pangenesis occurs: it is negligible for a flat direction carrying L v 
but no B v , and is m p x q x vis ~ 1 GeV x q x vis for a flat direction carrying B v but no L v . 
The above considerations yield a (2.5 — 3) GeV bound, as denoted in Eq. (4.74). 

The inequality (4.74) is satisfied for 

where we used Eq. (4.72). As per Eq. (4.67), Q-balls with charge larger than 

Q > Qstable « 10 32 (f^) 2 ( 4 -76) 

are thus completely stable. We mark this regime in the plots of Fig. 3. 
Medium-sized Q-balls with mass-to-charge ratio in the range 

3 GeV< — ^ <m NLSP ~m s . (4.77) 

q x y 

As per Eqs. (4.72) and (4.67), this occurs for 

70< fe 6K< 10 4(_=|_), (4.78) 
and Q-ball charges in the interval 



(4.79) 



The decay of Q-balls in this regime can produce only i?-parity even particles and graviti- 
nos. (Some NLSPs may also be produced at the late stages of Q-ball decay, when their 
charge has decreased and the mass-to-charge ratio has exceeded ttt-nlsp-) This scenario has 
been studied in Ref. [90], where it was shown that the decay into gravitinos is suppressed, 
and that the gravitino abundance resulting from Q-ball decays is subdominant if the el- 
lipticity of the condensate is e > 10~ 6 (m s M m /10 12 GeV 2 ) ~ 10~ 3 (?n 3/2 /l GeV). For the 
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parameter A of Eq. (4.25), the ellipticity is given by Eq. (4.41), and always satisfies this 
bound. 

Thus, for the range of b given in Eq. (4.78), the only constraint from Q-ball decay we will 
consider is Tq > 10 MeV, as in Eq. (4.49). Using Eq. (4.73), this limit becomes 

5 



b VlGevJ < U MsOOGevJ ' [ * U) 



Small Q-balls with mass-to-charge ratio 

Mq 



Q 

As per Eqs. (4.72) and (4.67), this occurs for 



> m NLS p ~ m s ■ (4.81) 



b c M < 70 A- , (4.82) 
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or 

,, 1022 ( ^^y^. ( , 83) 



In this case, Q-balls can decay into .R-parity odd particles (besides the gravitino), in 
addition to particles that are iZ-parity even. This is similar to the b < 1 case, and we shall 
adopt the same constraint, Tq > m 3 / 2 /20, as in Eq. (4.64). Given Eq. (4.73), this becomes 



m( n^_y <2.10Vjf "° V . (4.84) 
VlGeW J Q\500 GeVJ V ' 



In summary, for the scenario of pangenesis explored in this paper, and for the GMSB 
case, we shall consider the parameter region bounded by the following relations: 



m s \ 1 



^ * 104 Isod^vJ V x ' (485) 

llGeVv 1 < ^UoOGevJ ' 1 j 



5 



6 



These boundaries arise respectively from: the opposite regime of Eq. (4.75), ensuring Q-ball 
decay; Eq. (4.80), ensuring Q-ball decay before BBN, at Tq > 10 MeV; and Eqs. (4.84) and 
(4.65), ensuring no overproduction of gravitinos from Q-ball decay. We sketch these regions 
in Fig. 3, together with the constraints on the reheating temperature given in (4.31). We 
emphasise that at least some of the above bounds can be relaxed, according to the preceding 
discussion. 
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Table 1. Charge assignments of dark sector fields. The partners A and A carry the opposite charges. 



5 Models of Pangenesis 

In this section, we present a minimal dark sector which allows for successful pangenesis and 
discuss its phenomenology. In subsections 5.1 and 5.2, we describe the particle content of 
the visible and dark sectors. We identify suitable FDs for pangenesis in subsection 5.3. In 
subsection 5.4, we discuss the stability of particles in association to the symmetries of the 
model. We describe the cascade of the anti-baryonic asymmetry to the lightest dark-sector 
particles in subsection 5.5. 

5.1 The visible sector 

We model the visible sector by the MSSM, which contains the chiral superfields Q, L, u c , 
d° and e c in the standard notation defined earlier. To ensure the cancellation of (B — L)- 
anomalies, we also include left-handed antineutrinos n c . At the renormalisable level, the 
symmetries then allow for the superpotential 



where, for simplicity, we have suppressed family indices and coupling constants, and omitted 
obvious SU(3) C and SU(2)^ index contractions. Note that dangerous £> v -violating terms are 
automatically absent due to (B — L)-invariance and that there is thus no need to impose an 
i?-parity. 

5.2 A simple dark sector 

The dark sector must contain at least one chiral multiplet A which is charged under B — L. 
To ensure anomaly cancellation and allow for a supersymmetric mass term, we also include 
a partner A with opposite charge assignments. In order for Eq. (2.2) to establish a relation 
between VM and DM, we have to ensure that the symmetric part of the DM can annihilate 
efficiently to massless (or very light) particles. This requires an annihilation cross-section 
somewhat larger than the weak scale [62] . The prospect of efficient annihilation of DM into 
SM states are thus severely limited [134, 135]. The interactions between the dark and visible 
sectors - due to the B—L vector multiplet and higher-dimensional operators - are constrained 
to be suppressed by a larger scale, and give considerably smaller annihilation cross-sections 
of the DM to the visible sector 15 . This motivates introducing an unbroken gauge symmetry 
in the dark sector under which the dark relic baryons are charged, allowing the symmetric 
part of DM to annihilate to the dark force carriers. 16 The resulting annihilation cross section 

15 If the gauge-interaction which establishes the connection between the visible and the dark baryonic asym- 
metries couples to the visible sector as B — xL with x < 1, the annihilation cross-section of DM into SM 
fermions could be significantly larger than in the case of B — L. The collider constraints on a gauged B — xL 
get significantly relaxed as x — > [136] . Such an interaction would of course require the existence of exotic 
fermions to cancel the anomalies. 

16 Alternatively, the DM could annihilate to light dark fermions via Yukawa couplings. 
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Table 2. Monomials up to 5th order in visible-sector fields which carry B — L but are otherwise 
singlets (taken from [95]). 

of the DM is a ~ a^m^M, where m DM is the DM mass and a D the fine-structure constant of 
the dark gauge force. As we shall discuss in Sec. 6 (see in particular Eq. (6.7)), cosmological 
constraints on our model require that a D > 0.1. Furthermore, the dark matter mass in 
pangenesis is predicted to be m DM = O(GeV) (see Eq. (5.35) in Sec. 5.5). Together, this 
guarantees a sufficiently large annihilation cross section of the DM. 

We shall consider the simplest choice for this gauge symmetry, an abelian U(l)z), under 
which A and A are charged. In order to allow for U(l)£)-invariant monomials with nonvan- 
ishing (B — L)-charge (which combine with visible-sector monomials to form suitable FDs as 
we discuss in the next section), we also introduce two chiral multiplets A and A with equal 
and opposite charges under U(1)d. For defmiteness, we shall choose the latter multiplets to 
carry vanishing (B — L)-charge, whereas we assign (B — L)-charges ±g DM to A and A, respec- 
tively. Our conclusions will not depend on this choice, though, as only the (B — L)-charges 
of the U(l)£j-invariant monomials AA and AA are relevant for the phenomenology of this 
model. The charge assignments are summarised in Table 1 and allow for the superpotential 

W d = ms AA + m x AA . (5.2) 

The masses ms and m\ are taken to be around the GeV scale. We outline how this mass 
scale may be generated dynamically in Sec. 5.4.2 and will present a more detailed model 
elsewhere [137]. 

The Weyl fermions in the supermultiplets A and A combine to form a Dirac fermion 5 
with mass ms- Similarly, another Dirac fermion A with mass m\ stems from A and A. We will 
denote the scalar superpartners of these fermions by the same symbols as the supermultiplets, 
i.e. A, A, A, A, and the U(l)£>-gaugino by \d. We collectively refer to the latter particles as 
the dark-sector superpartners. 

5.3 Combined flat directions 

We shall now discuss how the visible and the dark sector particle content of the previous 
sections allows for flat directions in the scalar potential along which pangenesis can be im- 
plemented. 
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Directions in field space are 'flat' if both the D-term and F-term contributions to the 
scalar potential vanish. Identifying FDs is simplified by the following correspondence [138]: 
Every gauge-invariant holomorphic monomial in the chiral superfields gives rise to a direction 
in field space along which the D-terms cancel. Let us consider a simple example: The D-teim 
which belongs to the \J(1)d vector multiplet reads 

D D = |A| 2 — |A| 2 — |A| 2 + |A| 2 . (5.3) 

Along the direction A = A = and A = A = ^> in field space, where the scalar <p parametrizes 
the direction, this D-term vanishes identically. Since A and A carry charge only under U(l)£>, 
all other D-terms vanish similarly if we set the remaining scalar fields to zero. This direction, 
which we associate with the monomial AA, is thus 'D-flat'. 

The F-terms lift some of the D-flat directions in the field space. Efficient asymmetry 
generation via the AD mechanism requires the F-term contribution to vanish at least at 
the renormalisable regime. 17 This imposes additional constraints in determining the flat 
directions of the scalar potential. For our simple example, the direction AA, this requirement 
is in fact not satisfied due to the mass term in the superpotential in Eq. (5.2). For the purpose 
of the AD mechanism, however, terms of positive mass dimension in the scalar potential are 
not problematic since they are suppressed by powers of the mass parameter over the field 
VEV in the early universe compared to quartic terms. (Such terms arise anyway from SUSY- 
breaking.) We will therefore still call such directions 'flat'. A given FD is typically lifted by 
higher-dimensional operators in the superpotential. The FD associated with AA, for example, 
is lifted by the operator Wm, ^ (AA) 2 /M,,. 

Let us now identify FDs which are suitable for pangenesis, for the model specified in the 
previous two sections. At the renormalisable level, B — L is separately conserved in the visible 
and the dark sector, producing the accidental symmetries (B — L) v and B^, respectively. 18 
Pangenesis is effective along FDs which carry charge under X = {B — L) v + B& but which 
are singlets under the orthogonal combination B — L (see Sec. 3). Since the chiral multiplets 
of our model are charged either under (B — L) v or under B&, such FDs necessarily involve 
fields from both sectors. 19 We denote by monomials of visible-sector fields which carry 
{B — L)-charge q but which are singlets with respect to all other gauge symmetries. Such 
monomials, up to 5th order in the fields, are listed in Table 2 [95]. Similarly, we denote 
corresponding monomials in dark-sector fields by of\ The set of such monomials can be 
straightforwardly constructed using Table 1. For g / and to lowest order, these are A A 
and AA. FDs which carry charge under X can be parameterised by products of these two 
types of monomials: 

Of) ■ . (5.4) 

Gauge-invariance of this product ensures the cancellation of D-terms along the associated 
direction in field space. 

The constraints from F-flatness have been systematically studied for the MSSM in 
Ref. [95]. Visible-sector monomials which fulfil these constraints at the renormalisable 
level are marked as 'F-flat' in Table 2. Note that Ref. [95] considered the MSSM without 

17 In fact, as recent progress in the AD dynamics has shown, for successful asymmetry generation, the 
directions should not be lifted by terms with dimension d ^ 4 in the superpotential (see Sec. 4). 

18 In the visible sector, B v and L v are in turn separately conserved, as noted earlier and will be discussed 
further below. 

19 The possibility of 'connector fields' which carry both charges was considered in [1]. 
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right-handed neutrinos. Including these fields yields additional F-term constraints which lift 
FDs involving the SU(2)i-invariant combination LH U already at the renormalisable level. 
Furthermore, the product of different visible-sector monomials which individually correspond 
to F-flat directions does not always correspond to an F-flat direction itself. This should be 
taken into account when using Table 2. Let us now consider monomials O^p . Since the 
dark-sector superpotential in Eq. (5.2) only contains mass terms, all such monomials fulfil 
the constraints from .F-flatness at the renormalisable level. 

In summary, FDs which are suitable for pangenesis are associated with monomials of 
the type shown in Eq. (5.4) with the constraint that the visible-sector component satisfies 
the .F-flatness conditions arising from the renormalisable superpotential of Eq. (5.1). The set 
of these monomials obviously depends on the (B — L)-charge of the 'basis monomials' AA 
and AA. For example, for the choices g DM = |,2, 3, such directions are, amongst others: 

(AA)Vtfcf for (?DM = i (5.5) 

AA {u c d c d c ) 2 for foM = 2, (5.6) 
AA d c d c d c LL for q DM = 3. (5.7) 

Since there are three families of quarks and leptons, each of these monomials actually gives 
rise to several FDs. The above monomials are (by construction) the lowest order gauge- 
invariant operators involving a particular set of fields. They ultimately lift the corresponding 
FDs, and they are responsible for the generation of a net X charge, according to the dynamics 
described in the previous section. 

However, it is possible that other non-renormalisable operators of lower order also con- 
tribute to the potential along the field directions associated with the monomials of pangenesis. 
This is the case if there exist gauge-invariant operators of the type 

WiiftD— ^xCn-iW , (5.8) 

where x is a field that does not participate in the FD under consideration, and O n -\((f>i) is 
a dimension-(n — 1) operator involving only FD fields fa (but not necessarily the entirety of 
them). To reproduce the dynamics described in Sec. 4, it is necessary that n > d, where d 
is the dimension of the A-violating monomial in the superpotential. This ensures that the 
term |(/)| 2 ( n-1 \ which the superpotential term of Eq. (5.8) contributes to the FD potential, 
has a negligible effect on the dynamics of the AD field. 

In the preceding discussion, we identified FDs for which no operators of the type of 
Eq. (5.8) with n ^ 3 exist. Let us now discuss operators of this type, potentially relevant for 
directions of pangenesis, with n > 3. There are no such operators involving only dark-sector 
fields, as long as the FD involves either the monomial AA or AA, but not both. Indeed, all 
dark-sector gauge-invariant monomials have the form 

where k and £ are positive integers. The contribution of these terms to the potential vanishes 
for k+l ^ 2 since A = A = 0, along a FD involving only AA (and vice versa if the FD contains 
AA). On the other hand, the richer field content of the visible sector may allow for operators 
of the type of Eq. (5.8), If such operators exist, they may reduce the multiplicity of the FDs 
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under consideration, or even lift them completely at a lower order than the monomials char- 
acterising these FDs. A systematic study of all non-renormalisable gauge-invariant operators 
is beyond the scope of this work. 20 We emphasise however that certain non-renormalisable 
visible-sector operators allowed by gauge invariance must be highly suppressed in order to 
satisfy bounds on the proton lifetime (we discuss this further in Sec. 5.4.1). In the following, 
we will thus ignore non-renormalizable visible-sector operators altogether, as their existence 
depends on unknown physics, and their absence can be accordingly justified. In this case, the 
monomials of the type in Eq. (5.4) are the leading terms in the superpotential which lift the 
associated FDs. These operators provide the required violation of X and CP for successful 
pangenesis. 

5.4 Stable particles 

5.4.1 Gravity- mediated SUSY breaking 

We will now determine the stable and metastable particles in the visible and the dark sector, 
beginning with the case of gravity-mediated SUSY breaking. All scalars and gauginos (and 
the gravitino) then obtain masses of order the soft scale m s . For simplicity, we shall focus 
our discussion on the case of integer (B — L)-charge assignments q DM for the dark-sector 
multiplets. It turns out that the set of stable particles then falls into two different classes, 
depending on whether q DM is even or odd. We shall first present the analysis for the case of 
even q DM and point out the differences for odd charge assignments at the end of this section. 

Let us first ignore interactions between the visible and the dark sector. Imposing {B—L)- 
invariance in the visible sector results in two U(l) symmetries, L v and B v , with the linear 
combination (B + L) v being an accidental symmetry (albeit violated by sphalerons). These 
symmetries stabilise, respectively, the lightest neutrino and the proton. In addition, the 
visible sector has a gauged U(1)em which keeps the lightest electrically charged particle - the 
electron - from decaying. All other gauge symmetries are either broken or confined. Fermion 
number conservation stabilises an additional particle. In order to identify this particle, it is 
convenient to absorb a discrete subgroup of {B — L) v into the fermion number and define the 
i?-parity of a given particle as usual by 

pW = (_l)3(B-i)v+2 8] (51Q) 

where s is the spin. Note that this symmetry is not imposed but follows from (B — L) v - 
invariance and fermion number conservation. As usual, all SM particles carry even .R-parity 
and all superpartners are odd. The lightest superpartner is thus stable. We will refer to this 
particle as the visible-sector LSP. 

The particle content and the superpotential of the dark sector allow for two U(l) sym- 
metries, which we can identify with and U(l)u. The 'dark baryon number' l?d guarantees 
the stability of 5 since it is the lightest particle charged under this symmetry. In analogy to 
the visible sector, we will therefore call 5 the 'dark proton'. The gauge symmetry U{\)o in 
turn stabilises A to which we will correspondingly refer as the 'dark electron'. In absence of 
gauge interactions, the superpotential also has four abelian continuous i?-symmetries. These 
symmetries are broken to a single U(l) i?-symmetry by the coupling to gauginos (and the 
gravitino) which in turn is broken to a Z2 imparity by gaugino masses and A-terms. Choosing 

20 For the MSSM, the operators which lift a given FD have been systematically determined in [95]. However, 
since our model includes right-handed neutrinos and (B — L)-invariance (which changes the set of higher- 
dimensional operators allowed by the symmetries), their analysis is not directly applicable. 
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vanishing U(l)i?-charge for 5 and A, these particles are even under the .R-parity whereas the 
dark-sector superpartners are odd. This symmetry thus stabilises the lightest superpartner, 
which we shall call the dark-sector LSP. The P-parity of a given particle can be written as 

= (_ 1 )-3S d +D+ 2s j (5 n) 

where D is the charge under U(l)£>. We thus see that this symmetry is a subgroup of XJ(1)d 
and fermion number conservation (the factor (— l)~ 3Bd was included for later convenience but 
does not affect the .R-parity assignments since all dark-sector particles have, by assumption, 
even B — L and thus even 

Let us now reintroduce interactions between the visible and the dark sector. These break 
the two separately conserved fermion numbers (or equivalently P-parities) in each sector down 
to one universally conserved quantum number. Using the fact that B — L = (B — L) v — B^, 
we can combine the two symmetries into a universal .R-parity: 

pfe«0 = (_ 1 )3(B-i)+0+2 S (5 12) 

This symmetry stabilises the overall LSP, i.e. the lightest particle among the visible-sector 
and the dark-sector LSP. The heavier particle of the two, on the other hand, decays via the 
(B— L)-gaugino or the gravitino, which couple to both sectors. For example, if the overall LSP 
resides in the dark sector, the visible-sector LSP can decay via an off-shell (B — L)-gaugino: 

LSP V — > \*b-l + SM particles — > LSPd + dark-sector fermions + SM particles . (5.13) 

The decay rate is suppressed by the large mass Mb~l of the (B — L)-gaugino. In order not 
to endanger BBN, we have to ensure that the lifetime is less than 0.01 s (see e.g. [139]). This 
gives the constraint 

L B-L ~ 



^lOV^y^)' GeV, (5.14) 

were f s is a phase-space factor (e.g. f s ~ 10~ 3 for a three-body decay) and we have assumed 
that the gauge couplings of B — L and V(1)d are of order 1. Note that decays between the 
two sectors via the gravitino are strongly suppressed due to its small coupling strength and 
are thus less important than the channel via the (B — L)-gaugino. 

In addition, the superpotential contains higher-dimensional operators which couple the 
visible and dark sector. For example for q DM = 2, the symmetries of the theory allow for the 
following X-violating terms to appear in the superpotential: 

W d ^-AA (n c ) 2 + ^3 AA (n c ) 2 [AA + AA + H u H d ] 



+ ^AA(LH u ) 2 + ^ I AAQQQH d n c + ... (5.15) 



Here the ellipsis refers to additional terms of order-7 and higher orders in the superfields 
and we have assumed that all X-violating terms arise at the same scale M* as the operators 
which lift the FDs suitable for pangenesis. The couplings in Eq. (5.15) break the separately 
conserved symmetries £> v , L v and B^ down to the global B — L. Since the lightest state 
charged under B — L is the neutrino, we expect that the ordinary proton, as well as the dark 
proton become unstable. 
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We shall first determine the lifetime of the ordinary proton which results for q VM = 2 and 
derive a bound on the scale M* from the requirement that the proton is sufficiently stable. 
Since B v - and L v -violating operators for g DM > 2 involve higher-dimensional monomials in 
visible-sector fields (cf. Table 2) and are thus suppressed by larger powers of M*, this bound 
ensures sufficient stability of the proton in these cases as well. 

The combination of the first and the last term of Eq. (5.15) gives the lowest order decay 
mode of the proton to SM particles via intermediate dark-sector particles, for example via 
the decay chain: 

p — > A* + A* + MSSM particles — > neutrinos + mesons . (5.16) 

This is the most dangerous channel involving the dark sector. As we will see in Sec. 5.5, 
all dark-sector particles are heavier than the proton and can thus only appear as off-shell 
states in the decay. Counting propagators of heavy scalars and gauginos, the decay rate of 
the proton resulting from Eq. (5.15) can be estimated as 

ryy\ 33 

(5 - 17) 

where m p is the proton mass and / p < 1 takes into account further suppression due to 
phase-space/loop factors and coupling constants. Requiring that the lifetime of the proton 
is larger than the observational bound, r v > 10 34 years, leads to a lower bound on the scale 
M*: 

M,>4GeV-^(^X) 11/5 . (5.18) 

This is always satisfied for the values of M* required to produce enough asymmetry. The 
bound is so mild because the decay proceeds via a dimension-8 and a dimension-5 operator. 

Note that the operators QQQL and u c u c d c e c , which allow the proton to decay involving 
only the visible sector, are allowed by the symmetries. If these operators would arise at the 
scale M* similar to those in Eq. (5.15), the proton decay rate would be too large, for the 
values of M* we considered in Sec. 4. In fact, these operators must be suppressed by a scale 
considerably larger than the Planck scale to allow for a sufficiently stable proton. A way to 
naturally obtain this suppression was presented e.g. in [102]. Since this issue already arises in 
the MSSM alone and is independent of our dark sector, we will not consider these operators 
further. 

Similarly, the higher-dimensional operators make the dark proton 5 unstable. The 
decay products include a dark anti-electron to ensure conservation of D-charge, and SM 
particles. 21 In Sec. 5.5, we will identify the dark proton with a component of the DM state. 
Accordingly, its lifetime should exceed the age of the universe, t v ~ 10 10 yr. Depending 
on the decay products however, there may be more stringent bounds on the lifetime coming 
from measurements of the cosmic radiation backgrounds to which the decays contribute. We 
shall now derive a bound on M* from the requirement that the dark proton is sufficiently 
long-lived if q DM = 2. In analogy to the case of the proton, this bound again guarantees 
sufficient stability for q DM > 2 as well. 

For q DM = 2, the first term in Eq. (5.15) allows the decay 

5^\ + 2v R , (5.19) 



21 This is assuming that mj > m\. In the opposite case, the dark anti-electron A would become unstable. 
For very degenerate masses m$ and m\, finally, both particles would be stable. 
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at 1-loop via sneutrino and {B — L)-gaugino propagators. The decay rate is 

r ~ /d 16^ M^V ' (5 - 20) 

* s B—L 

where fd ~ 10~ 3 contains the phase-space suppression of the 3-body final state. Because 
the neutrino states produced in this decay are sterile, no bounds from radiation backgrounds 
apply. Requiring > 10 10 yr translates to a lower bound on M*: 

^-^f^'WH^)- (6 - 2i) 

This is significantly stronger than the constraint of Eq. (5.18), but is still always satisfied for 
the values of M* required for successful pangenesis (see Figs. 1 and 3). 

The dark-proton decay mode of Eq. (5.19) is the channel with the largest decay rate 
since all other l?d- v iolating operators are suppressed by larger powers of M*. These higher- 
order operators induce DM decay into detectable SM particles (e.g. photons, charged leptons, 
active neutrinos), in which case constraints from cosmic radiation backgrounds apply. Even 
though these constraints on the lifetime of the dark proton are stronger than that imposed to 
derive Eq. (5.21), it turns out that the resulting bounds on M* are weaker than Eq. (5.21). 

In the discussion so far, we have assumed that q DM is even. Let us now consider the 
case of odd (B — L)-charge assignments in the dark sector. The dark proton 5 then becomes 
an odd state with respect to the universal -R-parity of Eq. (5.12). As the lightest imparity 
odd particle, it is absolutely stable. The dark electron A remains stabilised by its charge 
under U(l)z). However, the lightest superpartner (among both visible-sector and dark-sector 
particles) - which for even q DM was the lightest .R-parity odd state and thus stable - becomes 
unstable due to the higher-dimensional operators that break the global B v , L v and B^ down 
to B — L. The assignment q DM = 1 is not suitable for pangenesis, as it generically allows for 
the superpotential term AA n° and breaks the X symmetry at the renormalisable level. We 
shall thus consider the case q DM = 3. 

For q DM = 3, the terms allowed in the superpotential to the two lowest orders are 

W D AA (n c ) 3 + — !-j AA d c d c d c LL . (5.22) 

These operators induce the decay of the lightest superpartner. For example, the first term 
allows the decay of the L>-gaugino 

Ad -> {5 + A*) or (A + A*) 5 + \ + 3v R , (5.23) 

where the off-shell A or A decays at 1-loop via sneutrino and (B — L)-gaugino propagators. 
The decay rate is estimated to be 

4 2 7 

r ~ fr WW~ ' (5 ' 24) 

* B — L 



resulting in a lightest superpartner with lifetime 
r~5.10 16 yr(- 



M* \ 4 /0.1\ 2 /500 GeV\ 7 / M BL /g 2 BL \ 2 / iq- 



10 13 GeV J \g D J \ m s J \ 1 TeV 
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That is, for all the values of M* of interest for successful pangenesis (see Figs. 1 and 3), 
the lifetime is considerably larger than the age of the universe, t v ~ 10 10 yr. In bee. o.o, 
we shall see under what conditions the lightest superpartner is a subdominant component 
of DM. Note that for the decay mode of Eq. (5.23), no constraints from cosmic radiation 
backgrounds apply - even if the lightest-superpartner abundance is comparable to that of 
DM - since the only SM particles produced are sterile neutrino states. Higher-order operators 
produce decay modes with detectable SM final states, but the corresponding decay rates are 
suppressed by additional powers of m s /M* and no constraints relevant for pangenesis arise. 

The decay of the ordinary proton into visible-sector particles can proceed via the com- 
bination of dimension-8 and dimension-6 operators arising from the terms of Eq. (5.22). The 
resulting decay rate is then suppressed with respect to the q DM = 2 case at least by a factor 
of (rrip/M*) 2 , and the lower bound on M* is even more relaxed than Eq. (5.18). For odd 
charges q DM > 3, the bounds are further relaxed. 

5.4.2 Gauge-mediated SUSY breaking 

Let us now consider the case of gauge-mediated SUSY breaking. The mass spectrum depends 
on the set of gauge groups under which the messengers are charged. In the simplest case, 
the messengers couple to all gauge groups in the visible and the dark sector. As usual, the 
gauginos then obtain masses of order the soft scale m s at 1-loop and the scalars at 2-loop. 
The same particles as in the gravity-mediated case are stable with the exception of the LSP 
which is now the gravitino. The bound on the {B — L)-breaking scale M B _ L in Eq. (5.14), 
that we imposed to ensure that NLSPs decay sufficiently fast to LSPs, accordingly no longer 
applies. It is replaced by the bound on the reheating temperature for large gravitino masses 
that we have discussed in Sec. 4.2.2. 

Another possibility is that the messengers are charged under the SM gauge groups and 
~U(1) b-l but not under \J(1)d- The visible-sector and B — L gauginos again obtain masses 
at 1-loop and all scalars at 2-loop (the dark-sector scalars via (B — L)-gauginos). Naively, 
one expects that the mass of the U(l)£)-gaugino Xd is generated via 3-loop diagrams (with 
intermediate (B — L)-gauginos and dark-sector chiral multiplets). It was shown in [140], 
however, that the mass in this situation actually only arises at 5-loop order and can be 
estimated as 



Here g D and g B _ L are, respectively, the gauge couplings of U(1)d and \J(1)b-l- It is thus 
possible that Xd replaces the gravitino as the LSP. This may ease the constraints on the 
reheating temperature from gravitino overclosure and NLSP decays that we have discussed 
in Sec. 4.2.2. Since the coupling of Xd to dark-sector particles necessarily involves heavy 
scalars, it decouples at temperatures below the soft scale and is thus a hot thermal relic. 
To avoid overclosure of the universe, m\ D must therefore be smaller than 100 eV which is 
readily achieved. 

Finally, let us briefly mention the case that the messengers couple to the SM gauge 
groups but to neither \](1)b-l nor U(1)d- Scalar mass-squareds m\ s in the dark sector are 
then generated via intermediate SM fields and {B — L)-gauginos at 4-loop, leading to 




(5.26) 




(5.27) 
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This case is thus promising for the dynamical generation of the GeV mass scale for the dark 
matter in pangenesis. We will present such a model, where this mass is induced via the 
dark-sector soft scale m^, elsewhere [137]. 

5.5 Dark atoms as dark matter 

Pangenesis generates an asymmetry in the dark baryon number B^, which is carried by the 
components of the A, A multiplets. Since the dark gauge symmetry U(l) d remains unbroken, 
a compensating asymmetry resides in the A, A multiplets. 

The scalar components of the chiral multiplets decouple, decay and transfer the asym- 
metry into the fermionic degrees of freedom, according to 

A —> 5 + \d or A^ 5 + X* B _ L -> S + LSP V + SM fermions (5.28) 
A —?■ X + \o or A^ X + X* B _ L -»• A + LSP V + SM fermions (5.29) 

and similarly for A, A. The decay modes on the left assume that the scalars are heavier than 
the XJ(1)d gaugino, m^, m\ > m\ D . If this is not the case, the dark-sector scalars can still 
decay via an off-shell Xb-l, as shown on the right. Of course, the decay of the dark-sector 
scalars presupposes that none of them is the overall LSP. This is essential for the success of 
our scenario in the generic case when the transmission of SUSY breaking to the dark sector 
is unsuppressed. We will return to this point later in this section. 

The dark baryonic asymmetry eventually resides in the dark protons 8, which combine 
with the asymmetric abundance of the dark electrons A to form D-neutral dark atoms. The 
latter constitute the DM of the universe in our scenario. After 5, X become non-relativistic, 
their symmetric thermal populations annihilate away into Z' D bosons. We shall now calculate 
the relic DM mass density and compare it with that of VM, using the structural relation of 
pangenesis, Eq. (2.2). 

The conserved B^ and D charge-to-entropy ratios are 22 

r)[B d ] = q DM T) S = ffl , (5.30) 
V [D] =ns-rix = 0, (5.31) 

where q DM is the X or, equivalently, B& charge of the 5 fermions (and of the dark atoms). 
Equations (5.30) and (5.31) give 

1 v\x] 

Note that Eq. (5.32), and thus the relations that follow, would still hold for any i?d-charge 
assignment for the fields A, A which maintained qs + qx = 9dm- In fact, shifting the baryonic 
charges of A, A by equal and opposite amounts is equivalent to kinetic mixing of V(1)b-l 
with U(1)d- 

In the visible sector, the number density of protons over entropy density rj p is 

Vp = V [B V ] = a s V [{B - L) v ] = a s &, (5.33) 



22 We denote 77/ = [n(/) — n(f)]/s, for a field /. Since we are concerned with particle populations whose 
symmetric part is annihilated away, this is also the particle-number-to-entropy ratio for that field, 77/ = n(/)/s 
or ??(/) = -n(f)/s. 
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where the factor a s takes into account the possible processing of the {B — L) v asymmetry 
by electroweak sphalerons. If the asymmetry is released in the thermal bath before the 
electroweak phase transition, then a s ~ 0.35 [141]. However, if the reheating temperature is 
lower than the electroweak scale, and/or the condensate fragments into Q-balls that decay 
after the electroweak phase transition, then a s = 1. 

Using the above, the ratio of the dark and visible matter abundances is 

^dm _ ms \rjs\ + m\ \rj\\ _ (m s + m\)/q DM _ R f( , Q/n 

^VM m p \Vp\ 

Thus, the observed DM mass density is obtained if the mass of the dark atoms is 

m DM ~ 5 q UM a s m p ~ q DM ■ (1.6 - 5) GeV , (5.35) 

where m DM ~ mg + m\ (the binding energy of 5 and A is negligible). The DM mass range 
of Eq. (5.35) and its scaling with the baryonic charge of DM is a general prediction of 
baryon-symmetric models, quite independently of the detailed structure of the dark sector. 23 
Equation (5.35) also shows that as long as q DM > 1, the proton is kinematically disallowed 
to decay into the dark sector, as assumed earlier in this section. 

We now return to the point made earlier, that the decay of the dark-sector scalars is 
necessary for the success of our scenario: Assume that one of the scalars, e.g. A, is the 
overall LSP. Then, inverse decays of the type A + 6 — > decouple at T ~ m\ D /2 > toa/2, 
while A is still relativistic. At this point, A carries a significant portion of the baryonic 
asymmetry, shared via chemical equilibrium with the fermionic degrees of freedom: 77 a — T]g. 
After the decoupling of the gaugino inverse decays, the individual U(l)# symmetries of 
each chiral multiplet are restored (see Sec. 5.4). This means that 7/a will be conserved 
(independently of rjs), and will contribute to the relic DM abundance (the symmetric part 
of the A particles annihilates into massless Z' D bosons). If toa ~ m s ~ 10 2 GeV, Eq. (5.35) 
cannot be satisfied (unless q OM ~ 10 2 ), and the A relic abundance would overdose the 
universe. We note, however, that if the SUSY-breaking transmission to the dark sector is 
suppressed and produces soft masses of order ~ GeV (see Eq. (5.27)), a dark-scalar LSP 
can provide a viable scenario of pangenesis. 

In order to allow the dark atoms to constitute the dominant part of DM, we have to 
ensure that the LSP is underabundant. Its freeze-out abundance can be sufficiently small 
if the LSP resides in the visible sector and is mostly wino or mostly higgsino [142-145]. 
Alternatively, the LSP may be the U(1)d gaugino, which annihilates into the dark-sector 
fermions 5 or A via a scalar exchange of A or A, respectively. The cross-section for this process 
is o" ann « a 2 D /m 2 ds . Cosmological constraints set a D > 0.1 (see Eq. (6.7)), which provides for 
a larger annihilation cross-section than that of a visible WIMP LSP. Additional enhancement 
of the LSP annihilation arises if the dark-sector soft scale is suppressed compared to that in 
the visible sector, m^s -C m s (in which case the LSP naturally resides in the dark sector). 

6 Cosmological constraints on dark sectors 

We now discuss constraints on our dark sector arising from BBN, large scale structure for- 
mation and the Bullet cluster. 

23 If sphaleron-type effects exist in the dark sector, and process the Ba asymmetry, the DM mass prediction 
should be modified by an appropriate factor. 
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6.1 Dark radiation and BBN 

A key feature of asymmetric DM scenarios, including those of pangenesis, is the annihilation 
of the symmetric part of the dark matter plasma into dark radiation. This gives rise to a 
constraint from BBN, because the additional radiation from the dark sector will increase 
the expansion rate of the universe during this epoch and potentially endanger the successful 
prediction of the primordial light element abundances. 

It has become customary to quantify an additional radiative component of the thermal 
bath during BBN by an "effective number of extra neutrino species", 5N c s. Currently, 
data on the primordial abundances allows extra radiation equivalent to about one additional 
neutrino species: <5iV e g < 1 at 95% C.L. [146]. This additional radiation is consistent with 
WMAP determinations of the relativistic energy density at recombination, which stands at 
5N cS = 1.34±°Jg at 68% C.L. [147]. It thus suffices to impose the constraint SN cS < 1. 

The radiative energy density due to relativistic dark sector species is determined both 
by the number of such degrees of freedom, and by the temperature of the dark plasma. The 
latter is not expected to be same as that of the visible sector for temperatures below the 
decoupling temperature T^ ec , the lowest temperature at which the two sectors are in thermal 
equilibrium with each other. The histories of how entropy is released from massive species 
into radiation will, in general, be different in each sector. For our purposes, we wish the dark 
sector to acquire a temperature T d that is sufficiently lower than that of the visible sector T v 
in order to suppress the relativistic energy density contributed by dark radiation. 

The evolution of the ratio T d /T v is easily tracked through entropy conservation, which 
implies that 

9v{T v )T^ = g v (T dcc ) 

ffd(T d )T d 3 g d (T dcc ) ' 1 • > 

where g v (T) and g d (T) are the temperature-dependent degree-of- freedom counts in the visible 
and dark sectors, respectively, and T dcc is the common temperature at the time the two sectors 
decouple. Using this, the <5iV e ff < 1 limit becomes 

9A (T dcc )< 36 V4 (M|fT) . (6.2) 

The visible and the dark sectors are kept in equilibrium via the gauged B — L interaction. 
They decouple at temperatures around the B — L breaking scale, T dcc ~ M B _ L > Mew- 
Thus, assuming MSSM content for the visible sector, g v (T dec ) 213.5. For the dark sector 
of Sec. 5, g d (T dec ) = 20, and <7d,BBN = 2 at the onset of BBN (T v ~ 1 MeV), so the 5N e fi < 1 
requirement is fulfilled. 

6.2 Dark-matter recombination and self-interaction 

Our dark sector furnishes 'atomic' DM, meaning that the DM today is actually a hydrogen- 
like bound state of fundamental dark fermions due to the unbroken TJ(1)d gauge interaction 
which acts like an electromagnetic force in that sector. This kind of DM has now been well 
studied (see e.g. [21, 148]), and we simply quote previous results and explain how they apply 
to the present scenario. 

The first constraint arises from the need to have successful large scale structure (LSS) 
formation. To ensure an appropriately fast growth of structure in the universe, overdensities 
driven by gravitational attraction must develop in the dark sector sufficiently early. From 
general considerations, it is well known that structure will only begin to form when the 
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universe becomes matter-dominated. But density inhomogeneities also cannot develop before 
the recombination of the dark atoms, because the growth of gravitational inhomogeneities will 
be disrupted by the relatively strong U(1)d force. To ensure good LSS, it suffices, therefore, 
to require that the dark atoms recombine into D-neutral bound states no later than the point 
of matter-radiation equality. The criterion for when recombination can be considered to be 
effectively complete is that the ionisation fraction is at most 0.1. Requiring that the residual 
ionisation fraction after recombination is no more than this value leads to the constraint [21] 

while the necessity to have recombination complete before matter-radiation equality en- 
forces [21] 

M D >0.5keV^ , (6.4) 

where fi D is the reduced mass of the dark atom system, and a D is the U(l)c fine-structure 
constant. Both of these constraints are mild. 

Finally, observations on the Bullet cluster place an upper bound on the scattering cross- 
section a for dark atom collisions with themselves, estimated to be [149, 150] 

<lcm 2 /g. (6.5) 



The cross-section is related to the dark atom properties through [148] 

a ~ 4-irk 2 a 2 , (6.6) 

where 3 < k < 10, and a Q = (a D H D )~ l is the Bohr radius. Using this in Eq. (6.5) leads to 
the bound 

'k\ (I GeV\ /3.2 GeV\ * 
»d ) V 7T1 e 

which is the most stringent of the constraints derived in this section. 



(—)(— ) 



7 Signatures 

7.1 Collider signatures 

The generalised and gauged B — L interaction is the most generic signature of pangene- 
sis. Collider experiments currently constrain Z' to be heavier than about 500 GeV for 
a B _ L = g 2 _ L /47r ~ 10 -2 [136]. A massive Z' may thus be within the reach of the LHC. 
If produced, it is possible to measure its invisible decay [151-164]. If the latter cannot be 
accounted for by neutrinos, such a measurement will be significant evidence for a baryonic 
dark sector. 

7.2 Dark-matter direct detection 

Antibaryonic dark matter can be detected via Z' B 
tering cross section per nucleon can be as high as 

af_ L « 8 • ID"- cm 2 ^ 



exchange. The spin-independent scat- 
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where we used m DM = 5 GeV. Current direct-detection experiments are not sensitive to such 
low masses, and there are therefore no constraints from XENON100 [165] or CDMS II [166]. 
However, models of pangenesis can easily accommodate larger masses, if the baryonic charge 
of the lightest I?d-charged particle is larger (cf. Eq. (5.35)). We note that the direct-detection 
cross section may be higher if the gauge interaction which establishes the connection between 
the visible and the dark asymmetries couples to the visible sector as B — xL with x < 1 (see 
footnote 15). 

It is also possible that the dark \J(l)o force mixes kinetically with the photon. If it 
is broken, the DM state is the dark proton alone, and thus carries D charge. The kinetic 
mixing between the dark photon and the ordinary photon provides another direct-detection 
channel with cross section 



where we used m DM = 5 GeV. Such an interaction could account for the regions favoured by 
DAMA and CoGeNT [167, 168]. 

We note in passing that a broken U(l)u can support a successful dark-sector scenario 
for pangenesis if the Z 1 has sub-GeV mass, such that the DM particles can annihilate into 
dark photons. In order for the symmetric part of DM to be annihilated efficiently via U(l)x), 
the gauge coupling has to be g D > 0.05. The massive dark gauge bosons Z' D can decay into 
lepton pairs of the visible sector via the kinetic mixing to the photon. For Z' D mass about 
100 MeV, the kinetic mixing can be e ~ 10 -4 [169, 170]. This allows for a rapid decay of 
the Z' D bosons immediately after they decouple. No extra relativistic energy density is then 
present at BBN, and no long-range force among DM particles exists, which in turn means 
no astrophysical constraints. The kinetic mixing between the dark force and the photon can 
be probed in current fixed target experiments [169, 171]. 

8 Conclusions 

The similar present-day mass densities for ordinary and dark matter suggest that they have a 
common origin in the early universe, and that the DM density is due to a particle-antiparticle 
asymmetry that matches the situation in the visible sector. A common cosmological origin 
for both types of matter is a key feature of baryon-symmetric models, which see the universe 
as being symmetric under a generalised baryon number that counts both ordinary baryons 
and dark matter particles. In these scenarios, an antibaryon asymmetry that exactly cancels 
the visible baryon asymmetry 24 is hidden from view in the dark matter of the universe. This 
enforces a tight relationship between the number densities of VM and DM and helps to 
explain their observed similarity. The DM mass scale is then predicted to be in the few-GeV 
regime. 

In this paper we reported on a detailed analysis of one of the main ways of producing 
a baryon-symmetric universe: the mechanism of pangenesis, which uses Affleck-Dine dy- 
namics. We focussed on all relevant aspects of both gravity-mediated and gauge-mediated 
supersymmetry breaking: asymmetry generation, gravitino production, and Q-ball forma- 
tion and decay. We showed that viable regions of parameter space exist for both cases, and 
when considering gauge-mediation we left the messenger scale arbitrary. The main results 
are summarised in Figs. 1 and 3. 

24 Modulo sphaleron reprocessing. 




(7.2) 
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We also presented a class of viable dark sectors. Typical features include GeV-scale 
masses and a dark-sector analogue of electromagnetism that is required for annihilating away 
the symmetric part of the dark plasma. The dark matter today is then typically predicted 
to be 'atomic', meaning that stable DM particles form hydrogen-like bound states due to 
their Coulomb attraction. Constraints from big bang nucleosynthesis, large scale structure 
formation and the Bullet cluster were derived, leaving a viable region of parameter space. 

There are two important non-gravitational interactions between the VM and DM in 
pangenesis. The most characteristic feature is a Z' gauge boson that couples to both sectors. 
The simplest models see the Z' coupling to ordinary matter via B — L. This particle can be 
produced in accelerators such as the LHC, but it will have a substantial invisible width into 
dark-sector species, distinguishing it from standard Z 1 possibilities. The discovery of such 
a Z' would obviously be a spectacular advance in understanding DM and its putative rela- 
tionship to VM. Kinetic mixing between the ordinary photon and the dark photon provides 
another potentially important interaction between the sectors, and can be probed in direct 
DM detection experiments. Overall, the following would constitute evidence for pangenesis: 
supersymmetry, together with GeV-scale DM mass and a Z' having an invisible width that 
cannot be accounted for by neutrinos. 
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